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wave ordering, and metal-insulator transitions, of lead- and potassium-doped barium 



O ' bismuthate. We discuss some of the early theoretical studies of these systems. Our re- 
^ ■ cent theoretical work, on the negative-t/, extended-Hubbard model for these systems, is 



also described. Both the large- and intermediate-?/ regimes of this model are examined. 



• using mean-field and random-phase approximations, particularly with a view to fitting 
O . various experimental properties of these bismuthates. On the basis of our studies, we 
point out possibilities for exotic physics in these systems. We also emphasize the differ- 
;_j ■ ent consequences of electronic and phonon-mediated mechanisms for the negative U. We 
show that, for an electronic mechanism, the semiconducting phases of these bismuthates 
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We propose other experimental tests for this novel mechanism of charge transport and 
comment on the effects of disorder. 
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The materials BaPbi_a;Bia;03, and Bai_a;Ka;Bi03, derived from barium bismuthate, 
BaBiOs, have many remarkable properties |], |, |, H, |lOl- These include 

structural and metal-semiconductor transitions as x and the temperature T are varied; 
diamagnetic, semiconducting (CDW ordered?) phases even in the presence of a partially 
filled Bi 6s-0 2po- antibonding band, with optical and transport gaps (2 eV and 0.24 eV, 
respectively, at x = 0) differing by nearly an order of magnitude; and superconducting 
Tc's (13K and 34K, respectively, for Pb and K doping) that are 3 — 5 times higher than for 
other three-dimensional oxide systems with similar, low densities of states at the Fermi 
level. However, partly because of all the attention devoted to the high-T^ cuprates, these 
bismuthates have not been studied extensively, either experimentally or theoretically. 

We review here the experimental properties of lead- and potassium- doped barium bis- 
muthate, give a summary of some early theoretical studies of these systems, and then 
provide an overview of our recent theoretical work |10|, which uses the negative-t/, 
extended-Hubbard model for these systems. The principal goal of our review is to consol- 
idate what is known (experimentally and theoretically) about these remarkable materials. 
We have organised the review as follows: Section I is devoted to the phenomenology of 
lead- and potassium-doped barium bismuthate. Section II is devoted to a survey of the 
models that have been put forth for these systems and to some of the early theoretical 
work on them. Sections III and IV are devoted to our work on these systems based on 
the negative-f/, extended-Hubbard model: Both the laige-U (Sec. HI) and intermediate-?/ 
(Sec. IV) regimes of this model are examined, using mean-field and random-phase ap- 
proximations, particularly with a view to fitting various experimental properties of these 
bismuthates. Section IV ends with some concluding remarks in which we point out, on 
the basis of our studies, possibilities for exotic physics in these systems, the different 
consequences of electronic and phonon-mediated mechanisms for the negative U, and new 
experimental tests for them. 
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I. Phenomenology of Barium Bismuthates 

This Section contains a detailed account of the normal and superconducting properties 
of the barium bismuthates. We begin with an introduction (Sec.I.l), summarize the 
remarkable experimental properties of these systems (Sec. 1. 2), and then discuss various 
experimental observations and their interpretations (Sec. 1. 3). In Sec. 1. 4 we give a short 
account of the electronic structure and configuration of these systems as obtained from 
band-structure calculations. 



I.l Introduction 



Compounds containing elements of Group VA (e.g., Sb, Bi, etc.) have drawn consid- 



erable attention from chemists for a long time. According to their electronic config- 
urations {ns^p^ in their outer shell n; e.g., for Bi n is 6 and the electronic configuration 
is [Xe]4/^^5(i^°6s^6p'^, where [Xe] corresponds to filled shells of xenon) they should form 
compounds with oxidation states 3"'' or 5"*", as these would yield closed, outer-shell config- 
urations for the corresponding cations. Indeed, these elements generally form compounds 
with valences 3^ or 5^, skipping the intermediate valence 4"''. Thus chemists were in- 
trigued when they were able to prepare compounds in which these cations have a nominal 
valence 4"*". However, the actual valences of these cations (in these compounds) turned 
out to be controversial; e.g., the valence of Sb in the halides of the type A2^SbXg~, (X 
represents a halogen and A a monovalent cation like NH4): Originally these compounds 
were regarded as instances of tetravalent Sb. Finally, careful Sb-Mossbauer studies 
confirmed that the actual compounds are like A4^Sb'^+Sb^+X]^2^, i.e., the Sb^+ ions charge 
disproportionate into Sb^"*" and Sb^"*" ions. 

A similar controversy has existed regarding the valence of Bi in BaBiOs, which was 
thought to be an example of a compound with tetravalent Bi. Here also there is some 
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experimental evidence (see later) that Bi charge disproportionates into Bi^+ and Bi^+ 
configurations, thus skipping the intermediate valence 4+. (Since Bi is Mossbauer inactive, 
a direct verification, as in the case of Sb Halides, is not possible.) The renewed and 
vigorous interest in this system was motivated by the discovery that these systems 
show superconductivity when doped with lead, with a maximum transition temperature 
Tc of 13K. For a long time thereafter, this was the highest superconducting Tc observed 
amongst compounds containing no transition metal. More recently, it has been reported 
0, |, ^ |T3|, |T^, that, on doping with potassium (which replaces Ba), this system shows 
an even higher superconducting transition temperature (T^ ^ 30/^ at 40% K doping). 

There are, in fact, about fifteen elements (not restricted to Group VA) in the periodic 
table that show this extraordinary property of valence skipping [0, |18|. For example, Ga, 
In and Tl show valences of 1^ and 3^, Ge, Sn and Pb of 2+ and 4"*" and Bi and Sb of 3"*" and 
S"*". V and Nb form compounds with valence states 3"'" and 5^, Ti with 2"^ and 4"'". In all 
these cases the intermediate valence states are either unstable or metastable. The formal 
valence state of these cations determines, to a large extent, the shape and orientation 
of the other orbitals (principally of the anion) in the compound they form; hence any 
theoretical calculations and predictions of the chemical nature and physical properties of 
these compounds would necessarily require this information. In addition, the microscopic 
physics underlying valence skipping is extremely important for the construction of suitable 
theoretical models. 



1.2 Summary of Important Experimental Observations 



The compounds BaPbi-a^Bia^Os and Bai-a^K^^BiOs and those in which Pb is replaced 
by Sb [0], Ba by Sr, and K by Rb P, |, |13|, |l5l have many interesting properties 

@, 1, 0, [i,0,|i|,|oH3. 
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(1) BaBiOs is an insulator even though it has a half-filled Bi6s-02p(j band. 

(2) The alloys exhibit a metal-semiconductor transition as the concentration of the 

dopant, say Pb or K, is varied [|l], ||, ^, |[. 



(3) However, at all dopant concentrations these systems are diamagnetic ^ |22|, 

so their semiconducting phases are clearly not Mott insulators. 

(4) In their semiconducting phases these alloys exhibit two energy gaps differing by nearly 

an order of magnitude (2 eV and 0.24 eV in BaBiOa); the larger gap shows up clearly 
in measurements of the optical conductivity and the smaller one in the temperature 
dependence of the resistivity 0, ||, The transport gap seems to be inacces- 

sible to optical-refiectance, photoconductivity, or photoacoustic spectroscopy. 

(5) The semiconducting phases of the lead-doped alloys show many properties (e.g., two 

Bi-0 distances) that indicate the presence of charge-density-wave (CDW) ordering 

[|, I, g I, n, 111,13,11. 



(6) These alloys undergo various structural transitions as functions of the concentration 
X and the temperature T; the structures exhibited by these alloys - ranging from 
monoclinic to cubic - can be obtained by slight distortions of an underlying cubic 



perovskite structure [H, 26, 2£ . 



(7) Both the lead- and potassium-doped systems are type-II superconductors with T^'s 
that are 3 — 5 times higher than those for other three-dimensional oxide systems 
with similar densities of states at the Fermi level (calculated from experimentally 
obtained values of 7, the Sommerfeld parameter ||, ^ ^ ^). 

Thus these systems are quite unusual and a theoretical understanding of their prop- 
erties is a challenging task, further complicated by the limited amount of experimental 
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data available and by the lack of agreement among experimentalists on many of the ob- 
servations and their interpretations. 



Recently |Q a family of Bi-based layered superconducting compounds have been 
synthesized which are isostructural with the tetragonal cuprate superconductor La2Cu04. 
These compounds, Ba:2-x^xP^i-yBiyOA-y, show an insulator to metal transition as a 
function of potassium doping. The metallic phase shows superconductivity with a highest 
transition temperature of about 14K. We will not discuss these compounds in this review. 



1.3 Experimental Overview 

This Section contains an overview of the experimental data available for these bis- 
muthates. We have tried, as far as possible, to collect results which are backed by more 
than one group and which have been obtained by using single- crystal samples. 

1.3.1 Phase Diagram 

The structures of barium bismuthates can all be obtained from slight distortions of an 
underlying cubic perovskite structure shown in Fig.I.l(a). In this structure, Bi (or Pb) 
atoms occupy the corners of a cubic unit cell, O atoms the face centres, and Ba (or K ) 
atoms the body centres. There are octahedral clusters of O atoms around each Bi (or Pb ) 
atom. Slight rotations or breathing-mode distortions of these octahedra lead to a variety 
of noncubic structures. For example, the distortions that yield the pure BaBiOs structure 
are shown in Fig.I.lb. Here there are two distinct nearest-neighbour Bi-0 distances - 
2.28 A and 2.12 A: oxygen octahedra around one of the sublattices of the underlying 
cubic structure contract and those around the other sublattice expand. This can be 
interpreted in terms of a freezing of the breathing mode of the oxygen octahedra leading 
to a CDW: Tilts of the oxygen octahedra further lower the symmetry of BaBiOs. It is 
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not surprising, therefore, that both BaPbi-^BiajOa and Bai_^K2^Bi03 show rich phase 
diagrams with many structural transitions. Their phase diagrams in the temperature T 
and the concentration x plane are shown schematically in Figs. 1. 2a and 1. 2b, respectively 



H, HI 



In addition to these transitions both these bismuthates show metallic, superconduct- 
ing, and semiconducting phases (Figs. 1. 2a and 1. 2b), the transitions between which show 
up most clearly in transport properties. The orders of these transitions are not very clear: 
The metal-superconductor transition is most probably second order. Some workers have 
suggested that there is a first-order transition at finite temperature from the metallic 
phase to the semiconducting phase [|l^ (hence the hatched, two-phase region in the phase 
diagram of Fig.I.2a). From the available experimental data it is not clear whether the 
superconducting and semiconducting phases coexist along any phase boundary and, if so, 
what the order of the transition along this boundary is 0, P). However, the maximum 
value of the superconducting seems to occur at nearly the same value of x as the 
metal-semiconductor transition. 

Electron diffraction study of the structural transitions in BaPbi^^Bi^Oa has been per- 
formed on powder samples by Koyama and Ishimara [|l^ and by Minami fl^ recently. 
The phase diagram as a function of doping x that Koyama and Isimara get differs from 
the one reported by Cox and Sleight [^] earlier. These authors do not get any or- 
thorhombic phase in their electron diffrcation studies. They find evidence for tetragonal 
and monoclinic phases at low temperature and a high temperature (~500K) cubic phase 
(the transition occuring at x = 0.4). They have also investigated the existence of any 
incommensurate phase and concluded in the negative, in agreement with Cox and Sleight. 

In early studies of the lead-doped compound [|T| it was suggested that superconduc- 
tivity could only occur in the tetragonal phase; however, subsequent work seems to sug- 
gest that it can also be obtained in specimens with orthorhombic symmetry p|. In the 
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potassium-doped compound, superconductivity has been observed only in the cubic phase 
^ . In both the potassium- and lead-doped systems, the superconducting phase appears 



right after the semiconductor-metal transition, with increasing lead or potassium concen- 
tration. 

1.3.2 Normal-State Properties 

Much more experimental data are available for BaPbi-^BiajOs 0, 0, ^ ^ 0, 



p5| , p6| , p^ , pEj , |3^ , |40| , ^ than for Bai-j^K^^BiOs. Some properties of the metallic 

phase (0 < X < 0.35 ) of BaPbi-ajBi^^Os are given in Table I.l: Resistivities are quite high, 
nearly independent of temperature, and increase from 210 fiQcm ( at x = 0) to 540/if2cm ( 
at X = 0.24). Carrier concentrations, obtained from measurements of the Hall coefficient 
39[ range from 1.42 x 10^°cm~^ ( at x = 0) to 2 x lO^^cm"^ ( at x = 0.2 ); the 



effective mass at x = 0.2 is 0.3me (obtained from plasma-edge data ^7\, |^). The specific 
heat C depends linearly on T at low temperature, but the coefficient of proportionality 
is rather low [^, ^ (~ 0.6 mJ/mole at x = ). However, plots of C/T vs 
show some deviation from linearity (requiring fits of the form C = 7T + f3T^ + aT^, 
with 7 = 0.6 mJ/mole K^, (3 = 0.42 mJ/mole K*^ and a = 1.1 mJ/mole K^). In the 
range < x < 0.2 optical-refiectivity data can be fit to a Drude form as expected in 
a conventional metal. Plasma frequencies, as obtained from optical experiments ||43| , 



increase with increasing x (i.e., with increasing Bi concentration) in the above range; this 
is consistent with the increase in carrier concentration with increasing x (Table I.l), as 
obtained from Hall-effect measurements. LAPW, band-structure calculations ^ 
are fairly successful in accounting for the observed properties of the metallic phase in the 
range < x < 0.2 (Sec. 1. 4); and the most recent band-structure calculations have 
also obtained reasonable results for the insulating phase. 

As the metal-semiconductor transition is approached by increasing x, the metallic 
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phase starts showing unusual transport properties 0, ^7\, ^ ^ (see Sec. 1. 3. 6): the tem- 
perature coefficient of resistivity becomes negative around x = 0.3 (Fig.I.3a) even though, 
at low-enough temperatures, the system becomes superconducting for x < 0.35, and only 
beyond x = 0.35 is true semiconducting behaviour (i.e., an activated behaviour over a 
large temperature range with the low-temperature resistivity showing no sign of a super- 
conducting transition) observed down to the lowest temperatures measured (Fig. 1. 3a and 
1. 3b). Reflectivity data show a gradual departure from Drude fits in the metallic phase 
from X = 0.25 onwards as the metal-semiconductor boundary is approached (Fig. 1. 4). 

Various transport and optical properties in the semiconducting phase of BaPbi_2.Bi^03 
have been measured, but there is a dearth of such systematic data for the Bai„2.K^Bi03 



system, though some measurements have started coming out [0, For BaPbi^^Bi^^Os 
the semiconducting phase extends over a wide range: 0.35 < x < 1. The temper- 
ature dependence of the resistivity in this phase yields a transport-activation energy 
Ea which goes from at x = 0.35, to 0.2 eV at x = 0.6, and finally to 0.24 eV at 
X = 1 (Fig.1.5). By combining the resistivity and Hall-coefficient data, the temper- 
ature dependence of the carrier concentration in BaBiOs (i.e.,x = 1) is found to be 
n(T) = 1.1 X 10'^'^[c'm~^)exp{—0.24eV/kBT). The pre-exponential factor is remarkably 
large and nearly equal to the number of unit cells per cm^, suggesting that the transport 
mechanism in this system is most probably intrinsic (and not extrinsic, i.e., not because 
of dopants, impurities, or defects). 

Recently some measurements of the resistivity of single crystals of Bai_a.Ki^.Bi03 have 
been carried out in the semiconducting phase ISO, ISTI] for different values of x. Hellman, et 



al. ffi their resistivity data for single-crystal and thin-ffim samples of Bai_a;Kij;Bi03 
to a variable-range hopping law p = poexp[(To/T)^/^], with a very high Tq (~ 3—6 x lO^K). 
According to them the states are very strongly localized and there is no indication of any 
nonlinear transport in these systems, so incommensurate CDW motion is possibly ruled 



9 



out. Resistivity data |B2[ for superconducting Bai_2;Kii,.Bi03 ceramics are also becoming 
available. These measurements seem to indicate that there is a variable-range-hopping 
regime with activation energy Tq ^ 10^ — 10^ in the normal state. It is now being argued 



by several groups ^ that, even in single-crystal samples of Bai^^^Ki^BiOs, two or 
more phases with different chemical compositions are present. (If so, then presumably 
one phase is metallic and the other semiconducting.) 

Sato et al., [^] have carried out a series of transport measurements on Bai_i.Ki2:Bi03 
thin films at x = 0.45 and found no evidence of contributions from impurities or inhomo- 
geneous distribution of potassium concentration in their annealed samples. The resistivity 
data for the best samples have a positive (but small) slope. Metallic temperature depen- 
dence is observed down from Tc upto 300K and the resisitivity at 300K is lower than 
500/if2cm. On the other hand negative temperature coefiiecient of resistivity was ob- 
served in samples that were grown at higher annealing temperature and these samples 
have resistivities much higher than 700/ificm. 

The Hall coefficient measurement on Bai_^Ki2,Bi03 thin films indicates that the 
charge carriers are electrons. This is in contradiction to the expectation that the carriers 
should be holes as monovalent potassium replaces divalent barium in Bai-^Kij^BiOs. The 
value of Hall coeffiecient Kh, (— )3 x 10^ cm^/C is about two times smaller than that 
expected from band structure calculations. The temperature dependence of Rh is quite 
small. These workers have also done Raman measurements and found that the height 
of the strong peak at 570 cm~^ that one observes at x = 0, is reduced and it broadens 
with potassium doping. This peak disappears close to the metal-insulator transition. The 
reduction of the peak as well as its broadening with doping is considered by these workers 
to be an indication that the remnant CDW like ordering that one finds in the metallic 
phase of the lead doped system Q (close to the metal-insulator transition) is absent 
in the potassium doped system. 
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Blanton, et al. [H^ have measured the reflectivity spectrum of polycrystalhne samples 
of Bai-^Ki^^BiOs as a function of doping x from x = to x = 0.4 in the frequency range 
from ~ 250 to 25000 cm^^. Their data on optical conductivity (obtained from the 
reflectivity data) indicate that the CDW conductivity peak at 16000 cm^^ at x = shifts 
downwards in energy and broadens substancially with doping. The spectrum at x = 0.4 
shows a broad contribution to the conductivity that is centred around 4000 cm~^. Blanton, 
et al. claim that this broad peak, has evolved from excitations across the single-particle 
(CDW) gap of BaBiOs and is an indication that some remnant CDW order persists even 
in the metallic phase close to the transition as in the case of lead-doped barium bismuth 
oxide Bozovic, et al. |^3[ have measured the far- infrared reflectance spectra in 



thin film samples of the superconducting Bai-^^Kia-BiOa at x = 0.4. They point out some 
overall similarity between the reflectance spectra of the cuprate superconductors and the 
bismuthates at optimal doping (corresponding to the highest Tc). All these spectra are 
nearly linear in frequency R{u!) ^ A — Buj (where A, B are constants, A ^ 1, B/h ^ 1 
eV~^) below the plasma edge. 

1.3.3 Charge Disproportionation and CDW ordering 

The actual valence state of Bi in the barium bismuthates is a matter of some dispute 
(Sec.I.l). Evidence in favour or against valence skipping by Bi in these systems is indirect. 
In particular, there has been some controversy in the experimental literature concerning 
the presence of a CDW in the semiconducting phase of these compounds 0, ||, |^, |2^ . 



The presence of the two valence states, Bi^+ and Bi^+, is intimately connected to this 
question. Many consider that the controversy has been settled more or less in favour of a 
CDW in the lead-doped system, but the results for the potassium-doped system are not 
clear yet. It is important to note that, to support a CDW in these systems, the charge 
contrast need not necessarily be 2e 0; indeed it is almost certainly less than 2e, as 
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suggested by recent, high-energy-spectroscopy experiments The CDW controversy 
in these systems is comphcated further by the methods used for sample preparation 0| : 
oxygen vacancies appear quite frequently in samples prepared in low-oxygen atmospheres, 
leading to inhomogeneities that make it difficult to identify long-range CDW ordering 
experimentally. Even in the lead-doped system, it is still not clear how the CDW order 
disappears with increasing inhomogeneity (or whether there is short-range CDW order 
away from perfect stoichiometry) . Optical-reflectivity spectra show that the CDW is 
quite robust (see below) and persists throughout the semiconducting regime (upto a lead 
concentration of 65 %). 

Neutron-diffraction (powder) studies p|, of BaKi_^Bia;03 show that the frozen, 



breathing-mode distortion of the oxygen octahedra leads to long-range CDW order in 
some of the semiconducting phases of BaKi.^^Bi^Os, but not in the semiconducting phase 
abutting the superconducting phase; so the cause of the semiconducting nature of this 
phase remains a mystery P, |10|, 

1.3.4 Optical Properties 



A number of optical-reflectivity measurements |0, |2^, ^ have been performed 
on the lead-doped system, both in the metallic and semiconducting phases, covering 
a wide range of frequencies, from far-infrared to visible. The metallic regime shows a 
characteristic Drude behaviour with a plasma edge, which shifts continuously towards 
high frequencies with increasing x. This edge is seen in the entire composition range 
without an abrupt change anywhere. Such a continuous variation of the reflectivity edge 
with X is indicative of the homogeneity of the system upto optical length scales. In 
the semiconducting phase the reflectivity spectra are dominated by low-frequency optical 
phonons. No discontinuous change is observed anywhere in the spectra in the metal- 
semiconductor transition regime. 
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In optical-conductivity spectra, obtained from reflectivity data, deviations from Drude 
behaviour occur as the transition is approached from the metallic side. This transition 



is marked by the transfer of the spectral weight into a high-frequency, quasi peak [g, ^ 



54| , pq] , which grows with increasing x into a peak containing most of the spectral weight 
in the semiconducting phase. The position of this peak shifts towards higher frequencies 
with increasing x as shown in Fig. 1. 6a. (Fig. 1. 5 shows this peak position as a function 
of x). It has been suggested that this peak (Fig.I.6a,b) has a profile which is like 
that of a CDW system with a gap along the entire Fermi surface (as in Ko.aMoOa) and 
not like the square-root absorption edge in conventional semiconductors. However, more 
recent work shows a broader peak The magnitude of the gap Eg, measured by the 

position of this peak in the optical- conductivity spectrum, is an order of magnitude higher 
than Ea, the transport activation energy P, ^ |lO|, e.g.. Eg = 2.0 eV and Ea = 0.24 
eV at X = 1, whereas in a conventional semiconductor Ea = Eg/2. The optical gap 
decreases with doping (Figs. 1.5 and 1.6) and, at the semiconductor-metal transition, it 
is quite small but nonzero. (The transport activation energy vanishes at the transition.) 
On the metallic side the plasmons show up strongly in the reflectivity spectra and mask 
the details of the optical gap. Throughout the semiconducting phase the optical gap is 
larger than twice the transport gap. The transport activation energy or gap Ea does not 
show up clearly in any optical, photoconductivity, or photoacoustic studies. 



Raman experiments p8|, ^ on BaBiOa show resonant enhancement and the ap- 
pearence of higher harmonics of the 569cm~^ mode phonon (the breathing mode described 
above), indicating that the optical gap Eg is related to the strong, electron-phonon in- 



teraction in the system ||28|, 0, ^ ^ 0, The resonant phonon energy for 



the potassium-doped system shows very similar behaviour and, as in the lead-doped case, 
the resonance decreases with increasing x. So far no one has reported any observation of 
two distinct peaks in the Raman intensity which could be identified with the two valence 
states of Bi directly. If there were a strong charge disproportionation in BaBiOs, one 
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would have observed a split peak with equal intensities. It is possible that the laser power 
used in these experiments creates oxygen vacancies and destroys the charge disproportion- 
ation locally. We have not come across any discussion of these issues in the experimental 
literature. 

1.3.5 Magnetic Properties 

Both the lead- and potassium-doped systems are diamagnetic in their normal phases 
over the entire temperature range of measurement and for all values of doping. The dia- 
magnetic contribution from the core, which constitutes almost the entire susceptibility 
in both the systems, when subtracted from the total susceptibility, leaves a tiny param- 
agnetic part which, in the lead-doped system, shows a somewhat unusual behaviour: it 
increases as the system becomes more insulating, over the entire range of x (Fig. 1. 7a). 
For the potassium-doped system the paramagnetic contribution remains nearly indepen- 
dent of X 1^, 1^, |2l|, |22|, The paramagnetic contribution from the valence electrons 
is very small and lies between ^ 0.8 to 2 xlO~^ emu/mol in barium potassium bismuth 
oxide. The temperature dependence of magnetic susceptibility is shown in Fig. 1. 7b. 

1.3.6 Metal-Semiconductor Transition 

The insulating nature of BaBiOa and the metal-semiconductor transition that occurs 
on doping with Pb or K cannot be understood on the basis of band-structure calculations 
0, ^ they yield a half-filled conduction band for BaBiOs. The simplest explanation 
for the insulating nature of BaBiOa invokes a CDW instability which opens a gap in the 
conduction band (see below). Both Pb and K doping should inhibit CDW ordering for two 
reasons: first, such doping moves the Fermi level down from its position at half filling (K 
doping removes electrons from the system; for Pb doping this follows from band-structure 
calculations); and it also introduces disorder. 
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We now summarize the experimental findings close to the metal-semiconductor tran- 
sition in the BaPbi_^Bia;03 system. There is a strong precursor effect of the metal- 
semiconductor transition as the Bi concentration is increased beyond x = 0.2. The 
resistivity data are shown in Figs. 1. 3a and 1. 3b. The temperature coefficient of re- 
sistivity becomes negative for x ~ 0.2. This would normally indicate that the system 
is already insulating. However, at low temperatures, the resistivity drops to zero upto 
X = 0.35 (Fig.I.3a), indicating a superconducting transition. Although the resistivity 
(for 0.25 ~ a; ~ 0.35) decreases with increasing temperature, it cannot be fitted to an 
activated behaviour over a wide temperature range. Such a fit is possible only beyond 
X = 0.35 (Fig. 1. 3b) where, down to the lowest temperatures measured, no superconduct- 
ing transition is observed. Hence the criteria used by experimentalists ^ for the metal- 
semiconductor transition in these systems are the absence of superconductivity at low 
temperatures and an activated behaviour of the resistivity over wide temperature ranges. 
There is then really no clear separation between the insulator and metal in the transition 
region and we might have a transition from a (Anderson-localized 1) semiconductor 
to a superconductor. 

Reflectivity spectra (Fig. 1. 4) show a marked departure from Drude behaviour from 
X = 0.2 onwards. The temperature dependence of the resistivity is unusual: From x = 0.2 
the resistivity becomes temperature independent and around x = 0.25 its temperature 
coefficient changes sign. The system is, however, metallic at that concentration (the 
resistivity does not show an activated behaviour). The thermopower varies linearly with 
temperature for x < 0.2; for x > 0.2 it no longer follows this metallic behaviour. 

Reflectivity spectra do not show any abrupt change across the metal-semiconductor 
transition in the high-frequency region. In the semiconducting phase there are several 
optic phonons in the low-frequency region. As the metal-semiconductor transition is ap- 
proached, the low-frequency conductivity is reduced and the spectral weight is transferred 



15 



to the high-frequency quasi peak as discussed earher. The quasi peak has a tail extending 
into the low-frequency regime indicating the possibility of states in the semiconducting 
energy gap. This tail disappears completely beyond x = 0.7. 



1.3.7 Superconducting Properties 

The superconducting phase occurs [|], ||, 0] in the concentration range 0.05 < 

X < 0.30 for BaPbi-^Bi^^Os and 0.37 < x < 0.60 for Bai-^K^BiOs. The highest is 
13K at X = 0.25 for the former and 34 K at x = 0.4 for the latter . The densities of states 



at the Fermi level for both these bismuthate systems are quite low |]T5|, ^ Given 
such meagre densities of states at the Fermi levels and the three-dimensional structures 
of these bismuthates, the superconducting transition temperatures are very high indeed. 
The specific-heat anomaly associated with the superconducting transition is very small 
Pm , Tunneling |^ ^ and infrared reflectivity measurements ^ yield 

2As/ ksTc ^ 3.5, the BCS value, where is the superconducting energy gap (measured 
from the chemical potential). For both these bismuthates there is a substantial isotope 
effect with oc M~", and a = 0.2 — 0.3 |2^ (A higher value of nearly 0.4 has been 



reported by Hinks, et al. 0, for the potassium-doped system.). The superconducting 



properties of these systems are summarised in Table 1.2. Note especially the coherence 



lengths ( ^GL ) which are nearly 70 — 80A, i.e., about 30 — 40 lattice spacings [T^, p4|. 

Sato et al. have measured the superconducting properties of Bai_2,.Ki2:Bi03 (at x = 
0.45) thin film annealed samples. Infrared reflectivity spectrum in the superconducting 
state at 5K (normalized by its value in the normal state at 30K) shows a clear peak 
at 60 cm~^ and a broad dip at 135 cm~^. The same features were earlier obtained 
by Schlesinger et al. in high density cintered samples with a transition width of 

4K. The superconducting energy gap that Sato et al. obtain is 3.7 ± 0.3 meV, in good 
agreement with Schlesinger et al. With the superconducting transition temperature of 
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their sample being 22.5 K, the value of 2A/A;bTc is about 3.8. Similar values have also 
been reported in the optical measurement point contact tunnelling measurement 
on Bai_2:Ki^Bi03 sintered samples at a; = 0.4 and tunnel junctions fabricated on 



thin films grown by molecular-beam-epitaxy and laser deposition method [^]. The 



conductance spectra for several tunnel junctions have been measured |^8|, ^ from 10 to 
100 mV bias voltage. The conductance spectrum in the region greater than 25 mV has 



no observable phonon structure p8| , |72| . This result seems to suggest that there is no 



strong electron-phonon coupling as has been assumed by many workers . 

The upper critical field in the barium potassium bismuthates has been measured in 
thin films of polycrystalline samples and in sintered samples [l70| , [7l[| . Unlike in the 
case of sintered samples, the thin films show a reasonably sharp resistive transition in a 
magnetic field. For the best thin film samples Hc2 shows almost linear temperature 
dependence down to 4.2K and seems to saturate aound 1.7K. 

The upper critical field (Hc2) measurements in polycrystalline samples of Bai-^^K^^BiOs 
by Crabtree ||6^, and Filippini et al., show that there is a significant upward curvature 
at fields beyond 10 Tesla. There is also noticeable flattening of the Hc2 curve close to the 
superconducting transition temperature. The value of the temperature coeffiecient of Hc2 
at the transition temperature {dHc2/dT at T^) for Bai_a;Ka;Bi03 presented in Table 1.2 
(~ —0.5 T/K) is for polycrystalline samples (from Batlogg, et al., p^)- Measurements 
on thin films |^8| indicate a much larger value of nearly -1.0 T/K (at x = 0.45). 

Hysteretic Josephson junction behavior has been observed in the I — V characteristics 
of multidomain films of Bai-j^Ki^-BiOa at temperatures upto 8K by Martin et al., |73| . 
They fit this behaviour to a model of Josephson junctions in series coupled by a shunt 
capacitor. 

1.3.8 The Effect of Oxygen-Deficiency on the Properties of BaPbi_a;Bi2.03 
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The effect of oxygen deficiency on the structural and transport properties of the lead- 
doped system has been studied by a few groups ||6^, |74|, |76|, Oxygen vacancies 
are easily created by sputtering in a thin film of a single- crystalline BaPbi_a;Bi03 sample. 
The principal effect of an increase in oxygen vacancies is to increase the resistivity in 
the metallic regime and decrease the superconducting Tc. Furthermore the resistivity 
of an oxygen- deficient metallic sample of BaPbi-xBixOs (at x = 0.3) goes down with 
increasing temperature and the sharpness of this fall in resistivity increases with increasing 
oxygen deficiency (Fig. 1. 8a). Thus, at a fixed temperature, as the oxygen deficiency 
increases, so does the resistivity. However, the "carrier density", as inferred naively 
from the Hall coefficient (for such a strongly disordered system), increases with oxygen 
deficiency (Fig.I.8b). The increase in carrier density with oxygen deficiency is thus not 



reflected in the resistivity results. Hence this suggests, that there is a strong carrier 
localization effect because of the defects created by the missing oxygen atoms, which more 
than offsets the increase in "carrier density" . 

It is surprising that the carrier density should increase with oxygen deficiency. For, 
as more and more oxygen is removed from the system, so are electrons, and one would 
expect a decrease in the carrier (electron) density (the Hall coefficient is negative). At 
very low temperatures (Fig. 1. 8b) it does seem that the more oxygen deficent samples have 
lower carrier densities; however, this trend is reversed at high temperatures. Enomoto, et 
al. argue that, at high temperatures, the carriers that were localized by the oxygen 
deficiency are delocalized (the activation energy for localization that fits their assumption 
is quite small ~ 18meV) and give rise to the enhanced carrier density; but the resistivity 
data, (Fig. 1. 8b) do not reflect this enhanced carrier density! This clearly calls into question 
the naive interpretation of the Hall-effect data. 



Another interesting observation |^3[ is that annealing a sample of the lead-doped 
system with x = 0.35 at 873 K and 1400 atm. (in an atmosphere of 80% Ar2 and 20% 
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O2) for 8 hours changes it from a semiconductor to a superconductor. This high-pressure 
anneahng decreases the resistivity markedly and, at 6K, there is a sharp superconducting 
transition with a width less than IK. However, if the sample is annealed at a pressure 
of 1 atm. (in an oxygen atmosphere) and at 600 K, it shows no superconductivity. This 
indicates that a metal-semiconductor transition can be induced in the lead-doped system 
by increasing the number of O vacancies. The sharpness of the superconducting transition 
is also greatly enhanced as the O deficiency decreases. A positron - annihilation study 
||77| on single- crystal samples of Bai_i,.Kj,.Bi03 in the range 15K to 300K shows that the 
sharpness of the superconducting transition and the value of the transition temperature 
depends strongly on oxygen vacancies, increasing as one approaches stoichiometry. There 
is no signature of positron localization and the positron lifetime is determined by the 
oxygen stoichiometry. 

A naive interpretation of these results is most probably fallacious. Firstly, there 



should be strong Anderson-localization effects in the electronic band, for, apart from the 
O vacancies, the concentration of Bi is about 30%. Secondly, the experimental system 
(at X = 0.3) is quite close to its metal-semiconductor transition and, as we saw earlier, 
this region is quite unusual. In particular, it is not clear whether the samples used in 
these experiments were metallic or semiconducting. Thus a more careful analysis of these 
results along with more controlled measurements are clearly necessary to understand the 
unusual observations reported here. 

1.4 The Electronic Structure of Barium Bismuthates 

The electronic structure of the barium bismuthates has been studied well using detailed 
band-structure calculations 0, 0, ^ as well as various experimental methods p8|. The 



electronic configuration of a Bi atom is [Xe]4/^'^5(i^°6s^6p'^ and that of a Pb atom is 
[Xe]4/^^5(i^°6s^6p^. Hence in the 4"*" valence state (which is the nominal valence state of 
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Bi and Pb in these bismuthates) Bi^+ has an open-shell configuration 6s^, whereas Pb'^^ 
has a closed outer-shell configuration. 

In both the lead- and potassium-doped systems, each Bi atom is surrounded by six O 
atoms which form an octahedron (Fig. I. la). There are ten orbitals in the "valence band" 



of the bismuthates [|^, ?5] given the BiOs unit cell, each O atom contributes three 2p 
orbitals, and each Bi atom contributes one 6s orbital (Fig. 1. 10). The Bi-0 bond lengths 
are the same for all the six O atoms in the cubic phase. We first examine the electronic 
structure of the cubic phase and later discuss the effects of the noncubic distortions. 

1.4.1 The Electronic Structure of BaPbi-^^BiajOs 

The electronic band structure of the bismuthates (Figs. 1. 11a - 1. 11c) has been com- 
puted by using the LAPW method The flat, low-lying bands are derived from 
Ba 5p states and the uppermost complex of unoccupied states are derived from Ba 5d 
and Bi 6p states. Note that the bands derived from Barium are admixed only marginally 
with the valence-band manifold in the middle. 

Of the three 2p orbitals at each O atom, one forms a a bond (Fig.I.lO) with the Bi/Pb 
6s orbital (spherically symmetric) and the rest form vr bonds. The vr — vr* manifold falls 
within the a — a* manifold (Figs. 1. 11 and 1.12). The Fermi level crosses the uppermost a 
antibonding subband (Fig.I.ll), leaving it exactly half fllled for BaBiOa (Fig. 1. 11c). As Bi 
is replaced by Pb, the flUing of this a* subband changes (Fig.I.lla-b). In addition, LAPW 
band-structure calculations show that Pb doping also introduces a noticeable chemical 
shift in the position of the 6s states. In all cases the Ba bands are well away from Ep. 
The principal features of these bands can be understood in terms of a tight-binding model 
P^ , |45| . Figure 1. 13a shows such a flt with a three-parameter, tight-binding calculation. 
The three parameters used are the orbital energies e^s = —4.1 eV, e2p = —1.9 eV and the 
nearest-neighbour, Bi-0 hopping parameter (along the a bond) tspa = 2.2 eV. As there is 
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no overlap in this model between Bi 6s and O 2p±, the entire vr — vr* manifold now shrinks 
to a set of perfectly flat and degenerate nonbonding bands. A five-parameter fit improves 
upon this, but the three-parameter fit clearly shows the tight-binding, a character of the 
important bands. We can, therefore, use this set of parameters to study the low-energy 
physics of these systems. This three-parameter, tight-binding model also predicts an 
electron-like, F-centered, simple-cubic Fermi surface for BaBiOs. The electron surface 
nests perfectly with an identical hole surface centered at the zone corner (R) (Fig. 1. 13b). 

Such a nested Fermi surface is susceptible to an electronic instability that, in the 
bismuthates, leads to a Peierls distortion, in which the O atoms surrounding every al- 
ternate Bi site are displaced towards them. This is the breathing-mode displacement of 
the oxygen octahedra, which doubles the unit cell and thereby opens up a gap along the 
entire, noninteracting Fermi surface, splitting the a* subband into a filled and an unfilled 
band. The effect of this Peierls distortion is shown, in a folded-zone scheme in Fig. 1. 14 for 
only the conduction-band states near Ep, using the three-parameter, tight-binding model. 
This gap is estimated to be nearly leV in a cubic BaBiOs structure, with alternate Bi-0 
bond lengths differing by 0.16A ||4^. This lowering of the energies of the states near Ep 
comes entirely from the Bi-0 bond alteration and hence stabilizes such a distortion. From 
the point of view of band-structure calculations, the role of doping is merely to inhibit 
the frozen-in commensurate distortion (since the Fermi level moves away from nesting), 
thereby producing a metal in which the conduction electrons are strongly coupled to the 
dynamic oxygen, breathing-mode phonons. 

Band-structure calculations work fairly well in the metallic regime of BaPbi-^Bij^Os 
(0 < a; < 0.2). The band crossing the Fermi level is derived from Bi/Pb6s-0 2p levels. 
Since Bi atoms supply one extra electron per atom to the conduction band, an increase in 
the Bi concentration pushes up the Fermi level and increases the number of carriers in the 
conduction band. The plasma edge shifts towards higher frequencies in this range exactly 
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as predicted by the band-structure calculations. The observed low density of states (DOS) 
at Ep is also brought out by calculations as the s-p character of the conduction band 
implies a low DOS at Ep (see below). The Fermi surface for x < 0.2 is nearly spherical 
as predicted by LAPW calculations. However, the band-structure results deviate from 
experimental observations from x > 0.2 onwards and cannot produce the semiconducting 
phase for x > 0.35. 

1.4.2 The Electronic Structure of Bai„a;K^Bi03 

Band-structure (LAPW) results are available for this system at various values 
of X. Figure 1. 15a shows these for BaBiOs using the observed lattice parameters of 
Bao.6Ko.4Bi03. The results are in close agreement with the ones obtained for the lead- 
doped system. The Fermi level crosses the uppermost of the ten-band, Bi 6s-0 2p complex. 
The flat, low-lying bands are Ba 5p states, whereas the uppermost complex of unoccupied 
states is derived from Ba 5d and Bi 6p states. The manifold of Bi6s-0 2p bands features 
a pair of broad (~ 16eV) cr and a* subbands extending from -12 eV to +4 eV at the 
zone corner. The a* band is half filled; so, as a K^"*" ion replaces a Ba^^ ion and removes 
an electron from the system, the Fermi level starts coming down the a* band, and the 
occupancy of this band is reduced to 1—x. Here also a three-parameter, tight-binding fit 
(similar to the Pb-doped case) works quite well. 

Recently Liechtenstein, et al. have carried out what is perhaps the best band- 



structure study of Bai-^K^BiOs. They use the LDA and the full linear muffin-tin-orbital 
method. For x = they find a lattice instability with a combination of both tilt {ut) and 
breathing {ub) distortions. This yields the experimentally observed monoclinic structure. 
For X = 0.5 they correctly find the cubic phase, with Ut = Ub = 0, to be stable. 

An interesting point that emerges from this band-structure study is that, when the 
details are put in, BaBiOs may have a rather small indirect band gap (Fig. 1. 15b). The 
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actual value of this gap, gleaned from Fig.I.lSb, is roughly 0.025 Ry or 0.34 eV; though this 
may not be very reliable. Hence one must consider the possibility that the small transport 
gap (0.24 eV) reflects this aspect of the band structure. However, optical experiments 
must be able to see this as a phonon-assisted threshold, starting around 0.68 eV, in order 
to be consistent with the transport gap. All the data that we have seen have thresholds 
starting at leV. (For recent spectroscopic data see Ref. |6^.) If the latter is taken as the 
indirect gap, the discrepancy with the transport gap is still too large to be understood in 
conventional terms. (For example, if one were to attribute it to polaronic or bipolaronic 
effects, as some groups have done ^ , the resulting Franck-Condon factors would be 

too large, leading to self-trapped rather than mobile carriers as required by the observed 
mobilities.) Thus one clearly needs an electronic negative-?/ mechanism of the sort we 
have discussed elsewhere 0] and which we will discuss below. (Of course, for a proper 
comparison with experimental data and a better estimate of model parameters it seems 
crucial to include the details of the band structure.) 

1.4.3 Comparison Between the Electronic Structures of Bismuthates and High- 
Tc Cuprates 

It is interesting to compare the electronic structures of the high-T^ cuprates and bis- 



muthates. There are some similarities H, 45] which become apparent when viewed against 



the electronic structures of conventional 'low-Tc' oxide superconductors, such as TiO. 
There are also very significant differences between the two classes of high-T^ oxide su- 
perconductors. Hence a comparison between their electronic structures and the resulting 
types of ordering (spin or charge) helps us in understanding the microscopic interactions 
that lead to the observed properties of these systems. 

Both these high-T^ superconductors are characterized by strongly hybridized bands 
(Cu3d-0 2p for the cuprates or Bi/Pb 6s-0 2p for the bismuthates), in which the Fermi 
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level falls within the uppermost a antibonding subband (Fig.I.12a-b), making it nearly 
half filled. In the cuprates this band is formed by the a bonding between Cu 3d(x^ — y^) 
and the 2p orbitals of the four surrounding oxygen atoms, and in the bismuthates, by the a 
bonding between the spherically symmetric and spatially extended Bi/Pb6s orbitals and 
the 2p orbitals of the six surrounding oxygen atoms. Note that, in low-Tc oxides (e.g., TiO, 
Lii+a;Ti2-a;04), the metal d states lie well above the 0-2p manifold and Ep falls within 
the TT — TT* manifold (Fig. 1. 12b). This has two effects: First, as the metal-3d orbitals are 
well above the 0-2p orbitals (e.g., in TiO this separation is nearly 5 eV),the hybridization 
between the metal d and oxygen p orbitals is quite weak, so the states near Ep for these 
low- Tc materials are predominantly of d- character (which is consistent with the observation 
that oxygen stoichiometry has little effect on their electronic properties). Second, the 
TT — TT* mainfold being triply degenerate (formed by the hybridization between the triply 
degenerate t2g orbitals, viz., dxy, dyz and d^^;, of the transition metal and the 0-2p orbitals 
of the neighbouring oxygen atoms), the DOS at Ep, and hence the Sommerfeld constant 
7, is much larger for them than for the high-Tc oxide superconductors, where the a* band 
is derived from the overlap of 0-2p orbitals with nondegenerate Cu3d— — or Bi 6s 
orbitals. The well-known Tc vs. 7 curve for various superconductors (Fig.I.16) illustrates 
this feature: The cuprates and bismuthates have a very low DOS at Ep compared to their 
superconducting Tc's. 

We have already discussed that a tight-binding fit for the important bands works quite 
well for the bismuthates. In the cuprates a similar, tight-binding fit also works well. For 
La2Cu04 tdpa = 1-8 eV and = Cp = —3.2eV; only two parametrs are required here as 
the 0-2p and Cu- 3d levels are nearly degenerate in the tight-binding fit. 

Like the bismuthates, the cuprates also have a nested Fermi surface {two dimensional 
as opposed to the three-dimensional one in the bismuthates) at zero doping. A system 
with a nested Fermi surface is susceptible to various instabilities like spin-density wave 
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(SDW), CDW, etc. The nominal valence of Bi in the bismuthates is 4"*" and that of Cu 
in the cuprates is 2+, which means that Bi and Cu ions have just one electron in their 
outer shells in these compounds. Thus one should expect these oxides to show fully devel- 
oped moments in their insulating states; however, these two systems show very different 
magnetic behaviour: The cuprates show a (somewhat reduced) moment whereas the bis- 
muthates are perfectly diamagnetic and there is no indication of a moment (e.g., in ESR 
measurements). In addition, the insulating phase of the cuprates is antiferromagnetically 
ordered, but in the bismuthates it exhibits CDW ordering. 

These differences suggest that the low-energy degrees of freedom in the cuprates are 
the spin degees of freedom, whereas those in the bismuthates are the charge degrees of 
freedom. Thus the dominant interaction in the cuprates is the onsite Coulomb repulsion 
{repulsive Udd) between two electrons (or holes) at a Cu site and not the long-range part of 
the Coulomb interaction, which favours CDW ordering. (An effective (antiferromagnetic) 
spin model can be obtained from a strongly correlated electron system via the Anderson 
superexchange mechanism, if the correlation is repulsive.) However, in the bismuthates, 
CDW formation and the consequent quenching of the spin degrees of freedom suggests that 
the repulsive, onsite part of the Coulomb interaction is suppressed relative to the repulsive, 
long range Coulomb interaction. The suppression can arise from a strong electron-phonon 
coupling between the Bi atoms and the breathing mode of the surrounding O atoms 



|78| , |79[| or from an electronic mechanism . There are reasonably convincing arguments 



(see below) for the onsite interaction in the bismuthates being effectively attractive and 
dominant relative to the nearest-neighbour. Coulomb repulsion. Thus, at the simplest 
level, the cuprates should be modelled by a Hubbard model with strong, onsite repulsion 
whereas the bismuthates should be modelled by using an extended Hubbard model with 

onsite attraction and nearest-neighbour repulsion. 

The experimental overview presented above suggests that the bismuthates are qual- 
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itatively different from both low-Tc, oxide superconductors and high-Tc, cuprate super- 
conductors. Their superconduting and unusual normal-state properties make these bis- 
muthates interesting, and an attempt to understand them a challenging task. So far there 
has been no fully successful theoretical treatment of these systems for all ranges of dop- 
ing. However, there have been quite a few attempts, with varying degrees of success, to 
understand some of the experiments mentioned above as we discuss below. 



II. Model Hamiltonians 

In this Section we review the various models and theoretical treatments that have been 
put forth to understand some of the experimental properties of the barium bismuthates, so 
as to put our work in perspective. These theories typically concentrate on some subset of 
the available experimental observations. Furthermore, most of these theories assume that 
the important, low-energy physics of these bismuthates can be understood on the basis 
of a negative-?/, extended Hubbard model with an attractive, onsite and a repulsive, 
nearest-neighbour. Coulomb interaction. An important issue, which has much bearing 
on a theoretical understanding of these systems, is the origin of the attractive, onsite 
interaction. Two points of view exist in the literature: 



(1) that the attraction arises because of a strong electron-phonon couphng ||78|, |79[; and 

(2) that the attractive interaction is electronic in origin ||18| . 



We discuss these two points of view below and the theories that use them. 

In Section II. 1 we motivate and describe the one-band, negative-?/, extended-Hubbard 
model for the bismuthates, and discuss the origin of the negative U. We also describe the 
(well-known) mapping of this model onto a pseudospin model in the strong-f/ limit. In 
Sec. II. 2 we discuss some earlier theoretical treatments of these models and analyze them 
in the context of experimental results. 
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II. 1 The One-Band, Negative-?/, Extended Hubbard Model For 
The Bismuthates 



At the BaPbOs end of the T — x phase diagram of BaPbi_a;Bia;03 (i.e., in the region 
X < 0.2) most of the properties of these bismuthates are conventionally regarded as being 
consistent with the simple band picture of Mattheiss and Hamman (Sec. 1. 4). A wide 
conduction band is formed because of the hybridization of Pb (or Bi) 6s orbitals with O- 
2po- bonding orbitals. An approximate, tight-binding characterization of this band may 
be made by choosing E^siBi/Pb) = —4.1 eV, = —1.9 eV, and tgpa = 2.2 eV. The 
band is partially filled as Pb^"^ contributes no conduction electrons to it (whereas O^^ 
contributes six p electrons per oxygen to the entire s-p band and two p electrons per 
oxygen to the bonding band). The superconducting transition temperature Tc is low; 
however, it increases with x in a manner that is consistent with the increase (with x) of 
both A^(0), the density of states at the Fermi level, and Aep, the electron-phonon coupling 
constant |56[] . 

For the regime x > 0.2, it is clear that such a simple band picture is inadequate because 
electron-electron and electron-phonon interaction effects are very important. It has been 
suggested |^ |79[ that some of the essential features of these interaction effects can be 
incorporated in an effective, semi-phenomenological, multiband, negative-?/, extended- 
Hubbard model where only the 0-2po- and the Bi/Pb 6s orbitals are kept explicitly, 
namely. 



<ij>(7 icr ia 

+v, (4c..)(p]..p,.')- 

Here cj^. creates Bi 6s electrons at the i sites of a cubic reference lattice (at this 
level of modelling, noncubic distortions are regarded as higher-order effects which can 
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be incorporated at a later stage), pj is the destruction operator for oxygen electrons (in 
the a orbital), tspa is the hopping matrix element for the transfer of an oxygen electron 
from the p„ orbital to a nearest-neighbour bismuth (or lead) atom, Ui and Vi account for 
onsite and nearest-neighbour Coulomb interactions, respectively, and is the chemical 
potential for oxygen electrons (with respect to the bismuth 6s level). 



It is easy to show |jTO|, |T8[ that, in the limit tspa ^ Cp, this reduces to a negative- 
U, extended-Hubbard model involving the Bi atoms only (i.e., involving electrons in the 
Wannier orbitals, centered on Bi atoms, of the antibonding Bi 6s-0 2pcr band only), given 
by the Hamiltonian 



n = -t J2 i4aCja + h.C.) - f//2 ^ fli^fli^^ + V ^i<Tnja' - fij^^i^ ^ ■'^) 

<ij>a ia <ij>cra' i 

where < ij > are nearest-neighbour pairs of sites, fii = J^a cjcrCia, is the chemical 
potential, i = |tpsp/ep, V = Viz\tspa/ ^p\^ and U = Ui. However, the vahdity of this 
reduction for realistic parameters and the relationship between t, V, and U and tspo-, Vi, 
and Ui clearly needs further study, for there are limits in which it must fail: in particular, 
it must fail in the limit when most of the Bi atoms are replaced by Pb atoms, so that 
BaPbOs turns out to be metallic as it is. Nevertheless, this effective, one-hand model, 
which has been proposed by a number of authors ^ |79|, ^ |4|, |85|, |86| , 



can indeed serve as a starting point, at least in the semiconducting regime, as it contains, 
the effects of the electron-electron and electron phonon interactions (at the simplest level 
as discussed below). 

In BaBiOa, l/A^I]jo- < fii^ > is constrained to be 1 (A^ is the number of lattice sites) 
corresponding to a nominal Bi valence of 4"*". In the Bai-j^K^^BiOs system, two important 
modifications are necessary: 

(1) Since K^"*" replaces Ba^"*", one electron is removed from the system and 

'^/NJ2ia < 'iT'itT >= (1 — x), corresponding to a nominal valence of (4 -|- x)"*". (2) There is 
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a random Coulomb potential because of K^"*" sites, which act as effective 1 impurities. 



In order to model |]79[ the BaPbi-^jBiajOs system, we must introduce quenched, ran- 
dom, site variables qj = (or 1), if site j is occupied by a Pb (or Bi) atom; i.e., qj = x, 
where the bar denotes an average over the quenched randomness. U and V in Eq.(II.2) 
are assumed to act only on Bi sites, which can be modelled by replacing in Eq.(II.2) by 
hiqi. In addition, we must introduce a random, site-energy term, eo I]icr'^io-(l ~ Qi), with 
eo = 2.5eV, corresponding to the difference between the local orbital energies of Pb and Bi. 
Note only Bi sites contribute electrons to the antibonding band, so l/NY,icr < ft-ia >= x 
(compare the case of doping above); however, Pb sites do not give rise to a Coulomb 
potential. 

The justification for the negative-f/ terms in Eqs. (II. 1) and (II. 2) is that it is the 
simplest phenomenological way of modelling the valence- skipping phenomenon, i.e., the 
stability of the Bi^+ (6s^) and Bi^'^ (6s°) states relative to the Bi'^+ (6s^) state, that many 
believe is characteristic of the chemistry of these bismuthates. As regards its origin, 
however, there are two points of view: the first phonon mediated and the other electronic. 

The approach in which the negative U arises from an electron-phonon mechanism is 



similar to Anderson's |87| suggestion (in a different context) that an attractive, electron- 
electron interaction can be generated if one couples lattice vibrations to electronic degrees 
of freedom strongly. Anderson had shown that the quenching of the magnetic moment 
in the A- 15 compounds and some amorphous semiconductors can be explained if there is 
such an attractive interaction at each site. This idea was utilized by several authors to 
understand the barium bismuthates, as outlined below. 



Rice and Sneddon and later Jurczeck and Rice |£9l used this idea and obtained 
the negative-t/ model by coupling the excess charge (Bi^+ being regarded as the neutral 
configuration) on a conventional positive-f/o Bi site to the breathing-mode, optic phonon 
of the surrounding oxygen octahedron. Consider the BaBiOa end of the composition 
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range: The energy required to displace each O atom is 



Hij = -Cx'^ - gx{qi -qj), 

where qi = (1 — rij) is the excess charge on a Bi site, (g^ — qj) is the charge difference 
between two neighbouring Bi sites, x is the (breathing mode) displacement away from 
the equilibrium position of the oxygen atom, g is the electron-phonon coupling constant, 
and C is the spring constant of the oxygen breathing mode. The minimization of Hij 
(with respect to x) yields x in terms q^ and qj and thus eliminates the oxygen degrees of 
freedom. The complete Hamiltonian that now describes the conduction electrons in the 
Bi 6s band, including the long-range. Coulomb interactions between Bi sites is: 



<ij>cr <ij> icr ij^j 

The effective, onsite interaction Uefj = Uq — zg'^/C is thus strongly suppressed 
compared to Uq, and can even be negative if g is large enough [z, the coordination 
number, is 6 in the three-dimensional cubic lattice considered here). Furthermore, 
Vnn = {e^/^rnn + g^ / C) IS the largest of the repulsive energies, where r„„ is the nearest- 
neighbour distance (lattice constant). 

If we ignore the longer-range, repulsive interactions present in the Hamiltonian (II. 3), 
we get an effective, one-band, negative-?/, extended-Hubbard model involving Bi atoms 
only (Eq. (II.2)). We can incorporate the effects of the long-range part of the repulsive 
interactions roughly by replacing e^/er„„ by e^a/er„„, where a is the Madelung constant. 
In any case V in Eq. (II. 2) is treated as a phenomenological constant, which must be 
obtained by fitting experimental data. 

It has been suggested recently |jl8| that there is an electronic mechanism by which 



this attractive, onsite interaction can arise in the barium bismuthates. This mechanism 
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is connected with the physics of valence skipping, electronic polarizability and screening 
processes. There are nearly 15 elements in the periodic table in Groups III and V which 
have a propensity for skipping a valence state (Ga, In and Tl favour valence states 1+ and 
3^ and not 2+; Bi and Sb favour 3^ and 5^ and skip 4^ (Sec. I)). The intermediate valence 
states are either unstable or metastable. This was well known in chemistry and attributed 
to the favoured valence states corresponding to closed-shell electronic configurations: The 
effective, intra-atomic, repulsion energy f/„ for the n*'^ charge-state of an ion, f/„ = En+i + 
En-i — '2En, where En is the total energy of the atom in the n^^ charge state, increases 
along a period for most elements as n increases, since the sizes of the orbitals decrease with 
increasing n. However, for the valence skippers, it is clear from Table II. 1 (compiled from 
measured or calculated ionization energies |]TB|) that Un has a minimum at the valence 
that is skipped. This reflects the larger correlation energy of the neighbouring, closed-shell 
configurations. All the U values, including those at the skipped valence state, are still 



positive, as is necessary for an isolated ion. Varma argues that this situation changes 
drastically when these ions are in a crystal. The value of U that should then be used 
is the effective, onsite interaction after the screening effect of the environment has been 
accounted for. Consider BaBiOs: The formal valence of Bi is 4+, so the configurations Bi^+ 
and Bi^"*" are strongly charged. Also, each Bi atom is surrounded by the highly covalent 
and easily distortable O octahedra, thus nonlinear screening of the atomic charges is 
important. A negative value of U (for Bi in the valence state 4"*" ) results if the reductions 
in the energies of 5"*' and 3^ configurations are sufficiently large compared to twice the 
reduction in the energy of the 4+ configuration. Since the 6s° (Bi^+) configuration is 
most strongly charged compared to the 4+ and 3+ configurations, it is screened most 
strongly by charge transfer from the O octahedra to the 6p states of Bi. Density-functional 
calculations of Hamman (quoted in Ref. |l^) give i?(6s°6p^) — -E(6s°) = —43.3 eV. As 
shown in Table II. 1, the electrostatic energy cost for the transfer of an electron is less 
than that obtained by linear interpolation between U^^ and t/3+ by nearly 15 eV. Hence 
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it is argued []TB| that nonlinear screening can bring the effective down to a negative 
value. 

Thus both electron-phonon and electron-electron interactions can lead to phenomeno- 
logical, effective Hamiltonians of the negative-?/ type. However, these two mechanisms are 
very different when considered in detail. In particular, if the electron-phonon mechanism 
is involved, the effective, onsite, attraction U is negative only on energy scales smaller 
than the frequencies of the breathing-mode phonons. If, on the other hand, an electronic 
mechanism is involved, U can be negative over typical electronic energy scales (of order 
eV). This difference has detectable experimental consequences (Sec. IV). 

Models (II.1)-(II.3) have strong correlation effects and many parameters, and are dif- 
ficult to treat: Hence many workers ^ have tried to see whether the large-f/, limit 
{U ^ zt) of model (II. 2), which is much simpler (see below), is adequate for describing 
the properties of the Bai-^^K^-BiOs system. In this limit we can project out the singly 
occupied (Bi^+) states and keep only the empty (|0 > i.e., Bi^+) and doubly occupied 
(|2 > i.e.,Bi'^"^) states at each site. The effective Hamiltonian in this reduced Hilbert 
space is easily derived via second-order, degenerate, perturbation theory (Appendix A) 
or canonical transformations, and can be written in terms of the pseudospin operators Si 
(that operate on the unoccupied and doubly occupied states \iO > and \i2 >, respectively, 
at the site i): 

St = cl^cl = \t2><i0l 

Sr = a^c,^ = \zO >< i2\, {11 A) 

SI = l/2{ni -l) = {\i2 >< i2\ - \iO >< i0\)/2. 

The effective Hamiltonian is 
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K// = J E (StS- + H.C.) + KY: S^S^ -BY.{2St + 1), (//.5) 

<ij> <ij> i 

where J = 2t'^/U, K = J + 2V, and the pseudomagnetic field B is related to the chemical 
potential fi via B = fi + U/2 — zV. Note that in the experimental system Bai-^K^BiOs, 
< Sf > is fixed at — a;/2; i.e., it corresponds to the fixed-magnetization ensemble. 

An alternative representation, that is physically more transparent than the pseudospin 
representation, uses the hard-core boson operators bj^ = Sl^ = cl^c]^ and bi = Sj^ = Ciic^. 
Clearly these represent real-space, onsite-paired 6s^ electrons, or Bi'^+ states, and have 
charge 2e. It is often useful to think in terms of these hard-core bosons, so we will use 
both pseudospin and hard-core-boson descriptions in our discussions below. The corre- 
spondences between the pseudospin model (II. 5) and the extended-Hubbard model (II. 2) 
are summarized in Table II. 2. We will discuss the viability of these models mainly for 
Bai-j^Kj^BiOs. To put our discussion in perspective, we review earlier theoretical treat- 
ments briefiy. 

II. 2 Brief Review of Earlier Theories of Barium Bismuthates 

II. 2.1 The Theory of Rice and Sneddon 

In one of the earliest attempts to understand the mechanism of pairing in the bis- 
muthates. Rice and Sneddon |7^, |79| developed the idea of the condensation of electrons, 
coupled strongly to the optic, breathing-mode phonons, into a crystal of real-space pairs 
(CDW order). These real-space pairs develop off-diagonal long-range order (ODLRO), if 
the band width is finite, and thence superconductivity. 

Rice and Sneddon worked with an effective, classical spin model. They expressed the 
Hamiltonian (II. 3), in the zero-bandwidth limit, in terms of a spin variable (where 
TTii = Hi — 1) which assumes the values 1, 0, and —1 at each Bi site; nti = +1 and rrii = —1 
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correspond, respectively, to the configurations with two electrons (Bi^+) and no electron 
(Bi^"*") at the site i, respectively, whereas rrii = corresponds to single occupancy (Bi^+) 
(the double degeneracy of this configuration is ignored). At the BaBiOs end, the Bi6s 
band is half filled, so the chemical potential is —U/2. The deviation of the chemical 
potential from its value at half filling (= —U/2) acts as a magnetic field in this effective 
spin model. Hence the effective, classical, spin Hamiltonian (in the zero-bandwidth limit) 



is just the S=l, Blume-Emery-Griffiths (BEG) model [B^ 



H = U,ff/2j2m^^ + V ^ m,m„ {II. 6) 

i <ij> 

where Ueff/2 = ^{Uo — zg'^/C), V = g'^/C + e^Q;/er„„ (clearly V > 0), and a is the 
Madelung constant (Sec.II.l). The intersite term is taken to be short ranged here and the 
effect of the longer-range part of the interaction is included via the Madelung constant. 
(Note the difference between the derivation of this classical spin model, obtained from 
model (II. 2) by setting the hopping t = and ignoring the the double degeneracy of the 
Bi^+ configuration (up and down spin states), from our earlier derivation in Sec.II.l of 
the quantum spin model in the limit of large-?/, from a negative-t/, extended-Hubbard 
model.) 

The Hamiltonian (II. 6) has been studied extensively |Q and its phase diagram is well 
known. For Uefj < zV (always satisfied if f/e// < 0), the ground state is antiferromagnetic, 
with nii = ±1 on the two sublattices (i.e., the CDW phase of model (II.3)). If Ueff > 
zV, the ground state has mj = at all sites (i.e., a nonordered phase with all Bi sites 
equivalent). 

Rice and Sneddon use mean-field theory to obtain the phase diagrams of model (II.3), 
with t = 0, and the BEG model (II. 6). Fig. II. 1 depicts their phase diagrams in the 
T — Ueff/2V plane. Though the mean-field phase boundaries obtained for these two 
models do not coincide, the topologies of the phase diagrams are the same. 
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At finite temperatures there are two types of excitations from the CDW ground state: 
two-particle excitations {nii —Trii); and single-particle excitations {rrii = ±1 — > 0), 
which break pairs. In the hypothetical limit Ueff/2zV — oo, a transition from the 
CDW to the nonordered phase can occur only via the two-particle excitations; this is 
a second-order transition from a crystal of paired electrons to a liquid of these pairs. 
As Ueff/2zV increases and passes through zero, the single-particle excitations become 
dominant and there is a first-order transition, with a discontinuous jump in the number 
of paired electrons (at T = this transition occurs at Ueff/2zV = 1). For Ueff/2zV > 1 
no phase boundary is encountered, as T increases; only the nature of the disordered phase 
changes from a liquid of paired electrons to a phase with no paired electrons. 

Doping by Pb'^^ is modelled as in dilute magnetic alloys since Pb*^^ behaves like a 



rrii = 0, nonmagnetic impurity. In this picture [^], the replacement of Bi by Pb removes 
an electronically active site from the system. Concentration fluctuations are assumed to 
be suppressed because of long-range Coulomb forces. At this level of modelling, the effect 
of doping by Pb is incorporated in the BEG Hamiltonian (II. 6) by the reduction of the 
local coordination number z. For Ueff < 0, both terms (site and bond)) in Eq. (II. 6) 
favour real-space pairing, and there is a second-order transition to the nonordered state. 
However, if Ueff > 0, then these terms compete. The reduction of the average cordination 



number z decreases the importance of the bond term, so it is argued |78| that dilution can 
drive the transition first order. Their phase diagram (Bethe-Peierls approximation) in the 
T/V—z plane is shown in Fig.II.2 for three different values of Ueff/2V = —5, 1.2, and 2.3. 

The only ordered phase of the classical model (II. 6) is the CDW phase. In order 
to have a superconducting phase, one must consider the quantum limit (with transverse 
spin fluctuations), so hopping (i.e., tij) must be included. With nonzero tij, it is not 
clear whether this model has an insulator-metal or an insulator-superconductor transition. 
Superconductivity may develop via a transition from real-space, local pairing (a 
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CDW) to momentum-space pairing (a BCS superconductor), aided by dilution (i.e., an 
increase of the effective, single-part cle band width); this scenario is described below: 

In the X —>■ 1 limit (the BaBiOa end), the CDW is the ground state, and two-particle 
excitations propagate only through high-energy intermediate states, resulting in a very 
narrow band. With dilution we get Pb^^ sites, which are rrii = 0, local, single-particle 
excitations over the CDW ground state. As the concentration of Pb sites increases, these 
excitations can move because of the hopping tij [= 2{z — z)t], so doping increases the 
band width for the motion of single-particle excitations. 

In the X ^ limit (the BaPbOs end), there are very few Bi atoms. They behave as 
impurities in BaPbOs. If f/e// > they prefer to be in the 4"*" state i.e., rrii = {z = 
zx ^ 0, so the bond term in Eq. (II. 6) can be neglected) and act as donors of electrons. 
As X and hence z increases, a metal-insulator transition becomes possible. If, however, 
Ueff < 0, the Bi sites act as negative-?/ centres and charge disproportionation into Bi^+ 
and Bi^+ is possible; i.e., an electron binds to the Bi^"*" site because of the short-range, 
onsite. Coulomb potential. This cannot occur unless |?7e//|/2 ^ 2{z — z)t. Hence Rice and 
Sneddon argue that []7S|, with decreasing z (i.e., increasing the Pb concentration) and at 
T = 0, there can be a transition from a CDW insulator to a metal only if |f/e//|/2 is not 
too large. If tij increases relative to the electron-phonon interaction (in the conventional 
BCS theory, tij is much stronger than the electron-phonon coupling), then the transition 
to an insulator from a metal occurs at a larger value of z, (i.e., a higher Bi concentration is 
necessary to make the system insulating). If t^j = 0, each electron is localized at a site and, 
for a system with N electrons, there are 2^ states. The spin degeneracy of the unpaired 
state is removed in the metallic state (when tij 7^ 0). As a result the entropic correction 
to the free energy goes as T^, in the metallic state, as opposed to T in the unpaired 
state (with tij = 0). Thus the metal-insulator phase boundary moves up perpendicular 
to the T = axis in the T/V — z phase diagram (Fig. II. 2). At lower temperatures. 
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then, the metal can become a supeconductor. This is the usual BCS (momentum-space) 
condensation of paired electrons into a superconducting state. Since the superconductivity 
here is mediated by optic phonons, the attractive interaction is retarded and all energy 
integrals (e.g., in the Tc equation) are cut off by c^o, the frequency of the breathing-mode 
optic phonon. 

II. 2. 2 The Theory of Jurczek and Rice 

The transition from real-space to momentum-space pairing as a function of doping 



ffql described above was studied further |79] in the limit of finite and large bandwidth 



to understand why the CDW phase in the lead-doped system is so stable (till a Pb 
concentration of 65%). The mean-field study of a strongly coupled electron-phonon system 
(with parameters taken from band-structure calculations ^) showed that the "local" 
CDW state is indeed stable upto large doping levels. (Band-structure calculations (Sees. I 
and II. 3.1) had predicted a strongly hybridized antibonding Bi/Pb 6s-0 2po- band crossing 
the Fermi surface as in a metal, which was clearly untenable since the lead-doped system 
shows a metal-insulator transition at a Pb concentration of 65%.) This theory begins with 
the Hamiltonian (II. 3) and an electron-phonon coupling that is strong enough to make 
Ueff = f/o — zg^/C negative, which can lead to both CDW and pairing instabilities. We 
outline this study below (first for the undoped and then for the doped case.) 

Undoped Case 

Here a; = 1 so the (uniform) CDW order parameter is: 

< Ui >=< ^c\^Ci„ >= 1 +6cexpiiQ.r^, 

where Q = 7r(l, 1, 1) and be measures the deviation of the valence at a Bi site from its 
nominal value (4+). Bi sites with < hi >= 1 + bc and < hi >= 1 — be are called A and B 
sites, respectively. With this choice for the order parameter and a conventional mean-field 

37 



decoupling the Hamiltonian becomes 



Ho = t Yl 4aCja + Ac^hi exp iQ.r^, (//.7) 

<ij>a i 

where Ac = Upbc/2 is the CDW gap parameter, Up = U^ff + 22;e^a/er„„, and a is the 
Madelung constant. The charge-disproportionation response- function is 

i?(6c) = l/iV^expiQ.r^- / dE nj{E]bc), {11.8) 
and, at half filling, the Fermi energy is given by 

1^1/NJ2 dEn,{E;K), (77.9) 
j ■'-^ 

where nj{E] h^, the local density of states at the site j with energy E, obviously depends 
on the value of which must be determined self consistently. If we multiply both sides 
of the self-consistency equation < fii >— 1 -|- bcexpiQ.r^^ by expiQ.r^- and sum over all 
sites, it becomes R{bc) — be- The local density of states (LDOS) that satisfies this self- 
consistency equation can be obtained , as usual, from the imaginary part of the Green 
function: 



1 1 

rij{E) =< j\S{E - Ho)\j >= - lim7m < j\- —\j > 

TT T?-*0 E — IT) — Ho 

= ^lim I mGj{z)\^=E-iri- (77.10) 
Finally the self-consistency equation becomes the usual mean-field gap equation: 
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where the sum over momenta is restricted to half the Brillouin zone, is the number of 
sites in the lattice, tk = ^ Z^k e^^ '^, and a are nearest-neighbour lattice vectors. 

Doped Case 

Jurczek and Rice accounted for dilution (i.e., Pb doping. Sec. II. 1) by replacing nj with 
Pjhj, where pj = 1, if site j is occupied by a Bi atom, and pj = 0, otherwise; the difference 
eo, (— 2.5 eV Q between the onsite energies of Pb and Bi, is an additional parameter. 
The following ansatz is used for the inhomogeneous, CDW order parameter: 

<nj >=pj{l + bcexpiQ.rj), {11-12) 
i.e., for a site with a Pb atom < nj >= 0. Mean-field decoupling yields the Hamiltonian 



Ho = t c^Cja + eoY^fijil - Pj) + Ac{x)Y,njPjexpiQ.rj, {11.13) 

<ij>a j j 

where Ac(a;) = \{Ueff + 22e^a/er„„)6c(x) and we have displayed the x dependence of 
various quantities explicitly, as the order parameter has to be obtained self consistently 
for each value of x. For BaBiOs there are two sets of inequivalent Bi sites, A and B. With 
the introduction of Pb, sites in each of these sets can have either a Bi or a Pb atom. 
Hence R{hc) (Eq. (II. 8)) becomes 



R{hc{x)) = 7--^ r[x / dE nXZWi{E; b,) + {1 - x) dE nXZrb{E; b,) 

ZXyL — X) J —oo J—oo 

dEn^,{E-b,)-{l~x) dEnWZWi{E-X)]. {ILIA) 

-oo J—oo 

The LDOS (Eq. (11.10)) must now be configurationally averaged. Even a coherent- 
potential-approximation (CPA) calculation of this averaged LDOS {ha-bi, etc.) is difficult 
since there are four different averages. Hence an orthogonalized-moment method has been 
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used [^. This method uses the continued-fraction representation of the one-particle 
Green function Gi{z). Figures II. 3a and II. 3b show the self-consistent LDOS for the Bi- 
A(B) sites at Pb concentrations of 20% and 60%, respectively. Note that, on doping with 
Pb: (1) there is no real gap, but only a pseudogap, in the DOS at the Fermi energy; and 
(2) there are very few Pb states in the pseudogaps of Fig.II.3. As the Pb concentration 
increases, the Bi subbands approach each other: the upper Bi subband is pushed down 
because of repulsion by the Pb band (the difference between the onsite energies of Pb and 
Bi is eg, — 2.5 eV) but the lower Bi band is hardly affected. 

It was found that the local CDW is not stable (i.e., be = 0) if the total coupling 
constant A(x) = ^\Ueff{x)\N{Ep, be = 0) is less than a critical value \c{x). {N{Ef, be = 0) 
is the DOS at the Fermi energy in the absence of the CDW; and this positive, dimen- 
sionless, total coupling constant is an effective measure of the onsite attraction relative to 
the bandwidth in the noninteracting limit.) Also, the local CDW in the BaPbi-^^Bi^-Oa 
system was shown to be stable for a wide range (0.35 < x < 1) with reasonable values 
of A. In this local CDW phase the DOS is nonzero in the pseudo gap (Fig.II.3), so, to 
account for the semiconducting nature of this phase in BaPb1_1.Bi2.O3, Jurczek and Rice 
argued that the states in the pseudo gap obtained here could well be localized because of 
diagonal disorder (i.e., Pb doping). 



An estimate of the electron-phonon coupling constant has also been made ||T9[: Uejf 
can be separated into an attractive part zg'^/C, from the electron-phonon interaction, 
and Uq, the repulsive, onsite electron-electron Coulomb interaction. These two parts can 
be used to define the electron-phonon and electron-electron coupling constants Xep and 
Aee, respectively (A's defined as above). When Uq exceeds Uqc = zg^/C the effective, 
onsite, interaction becomes repulsive and the ground state exhibits a spin-density wave, 
not a charge-density wave (if near-neighbour repulsion is negligible compared to the on- 
site term). Hence, in the CDW state. Age is maximum at Uq = Uoc- This maximum 
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value combined with typical values of Uo and the value obtained for Ac gives an estimated 
lower bound for Agp: it is found that Agp — 1 (for cq — 1.7 eV); i.e., according to this 
theory, the system is just in the strong electron-phonon-coupling regime. (This should 
be contrasted with the large values obtained in the bipolaron-superconductivity theories, 
discussed below.) The values of both Aep and eo compare reasonably well with experimen- 
tal observations: eo — 2.5 eV (atomic estimate given above) and the experimental value 
of Aep — 0.6 in the lead-doped system. These values have been used to obtain the 
superconducting Tc via the McMillan-Morel-Anderson Formula (McMillan, (1968; Morel 
& Anderson, 1962) and with Aep(&c = 0.65) = 0.7: = x^^-P(ltt62X,^) ) • Here 

< a; > is the average phonon frequency (this is normally taken to be the frequency of 
the phonon that couples to the electronic degree of freedom or the Debye temperature 
of the system), fi* is the Morel- Anderson Coulomb pseudopotential (Morel & Anderson, 
(1962). Jurczek and Rice obtained Tc > 7 K with < u >^ 150cm^^, the frequency of the 
optic-mode phonon, and 0.067 < /i* < 0.116. 

The main results of the above analysis is that the insulating, local CDW phase of 
the lead-doped bismuthate is quite robust against doping. There is a strong electron- 
phonon coupling in BaPbi_a,Bi^03 responsible for the attractive electron-electron inter- 
action that leads to the CDW and superconducting instabilities in this system; and the 
electron-phonon coupling constant is fairly large. We will briefly discuss the merits and 
shortcomings of this model below. 
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II. 2.3 CPA On An Attractive Hubbard Model: Semiconductor-Super- 
conductor Transition 



Yoshioka and Fukuyama |^ have considered a negative-?/, Hubbard model with 
onsite disorder for the bismuthates. They believe that, in BaPbi^^-BixOa, there is no 
semiconductor-metal transition but only a transition from a semiconductor to a super- 
conductor (at T = 0). Although such an assertion is not very unreasonable (Sec. 1.3), 
there is no firm experimental evidence for it. In their model then, as long as there is a 
gap in the DOS due to the difference in chemical potentials of Bi and Pb, the system 
remains insulating. As soon as this gap vanishes, the attractive onsite interaction leads 
to a pairing instability, which, in the doped case, leads to a superconducting transition 
(at half filling both the CDW and superconducting phases are degenerate, but doping 
favours superconductivity). We review their theory briefly below: 

In the negative-?/ Hubbard model they consider, disorder because of the substitution 
of Bi by Pb is modelled by a site-dependent chemical potential. Site energies and 
€3 are associated with the atoms A and B (see below) respectively. There is an onsite, 
attractive interaction mediated by a strong, electron-phonon interaction. They further 
assume that atom A has no conduction electron while B has two (see below). So at 
x=l the conduction band is full and at x= 0, it is empty . Now the connection between 
this picture of a binary alloy and the experimental system BaPbi_2.Bi^.03 is made in the 
following manner: 

The hybridization of Bi/Pb6s and 2p orbitals creates ten bands near the Fermi 
level (Sec. I, Fig. 1. 11 and 1.13). The Fermi level falls in the highest a* subband which is 
half-filled at x = 1. A unit-cell-doubling distortion creates a gap along the entire non- 
interacting nested Fermi surface (Sec. I) crossing this subband. Hence this subband splits; 
and the lower split subbands is filled at x = 1 and empty at a; = 0. Such a distortion 
corresponds to the frozen-in, breathing-mode phonon observed in this system. If the 
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difference between Bi and Pb site energies were much smaller than this gap, then for 
discussing the low-energy physics of this system it would be sufficient to concentrate on 
the lower split subband only. 

In the bismuthates, actually this is not the case; the Pb site energy is higher and the 
difference between atomic levels of Pb and Bi is ~ 2.5eV) Then the mapping of the real 
system onto the alloy is done by assigning the PbOa cluster to A atom and the BiOa 
cluster to B atom (and ignoring the existence of the CDW for other purposes). 

The resulting model is given by the Hamiltonian 

<ij>cr icT 

with €i — €a{b) for A (B) atoms. The randomness is treated in the CPA, whereas the 
attractive interaction is treated via standard, mean-field approximation, it is i.e., assumed 
that < -hi >= nA{B) for i belonging to A(B); n^ and n^ are determined self-consistently. 

In this approximation the mean-field Hamiltonian becomes 

77 — ^ ] (tij -\- €i(jSij^c^^CjQ-. 

<ij>a 

where Qo- = |(ei - C/ < fii-a >) 

To study the effect of disorder, we need e, the difference between the chemical potentials 
of the A and B site in presence of the mean-fields, i.e., 

£ = £(0) - U{nA - ub) 

where e^^^ = ^a — ^b- The local densities of states at A and B sites are given in the CPA 

by 
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where the self energy calculated in CPA, is E(^) — z — 1/F{z) — F{z)/A and 

F{z) = llNY,G{KiuJn). 

k 

Now ua and ub are determined self consistently via 



/Ep 
dE pa{b) {E) 
-oo 



and the lever rule (1 — x)nA + xub = x. Finally, assuming that there is a gap in the 
density of states in the semiconducting phase, we get (at T = 0) 



Ua = -wVe^ and = 1 - - — -w'^/e^. UF17) 
4 4 

where, w is half the noninteracting bandwidth, and thence 



which gives the condition that s^^^ should satisfy to produce a gap in the density of 
states. Figures 11. 4a and II. 4b show, respectively, the boundary between semiconducting 
and metaUic phases and the DOS. The total DOS p{E) (Fig.II.4b) is obtained from 
p{E) = (1 — x)pa{E) + xpb{E), where x is determined self consistently from Eqs. (11.17) 
and the lever rule. With the DOS thus obtained, the superconducting Tc is calculated 
from a Cooper instability calculation in the framework of the CPA . The expression for 
Tc is formally the same as for a pure system. 

U n ^1=^ ,^.tanh[{E-EF)/2Tc] 



l^-T dEp{E)'- 

Z J —wn 



-WD {E — Ep) 

while p{E) is determined in the manner described above. 
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Fig. II. 5 shows Tc versus x for different values of U . The cut off uj£, is some phonon 
frequency whose value is taken to be larger than the bandwidth in this calculation, but 
the phase diagram qualitatively the same. If the cut-off is reduced, decreases. Hence 
depending on the parameter values of the system {w, and U) the model system be- 
comes either a superconductor or a semiconductor as x is varied. Note that the phase is 
symmetric about x = 0.5. 

The results obtained here do not really conform to the experimental observations in 
the BaPbi_a;Bia;03 system. The insulating phase here has no CDW order as observed 
in the system. This is a very simple-minded model only explore the possibility of a 
semiconductor-superconductor transition but it may not be relevant for the BaPbi_j;Bij.03 
system. 

II. 2.4 Bai-^K^BiOs As A Doped Peierls-Insulator 

The models we have discussed so far are for the lead-doped bismuthate. They were all 
put forward before the discovery of superconductivity in Bai-^^K^-BiOs. Subsequent theo- 
retical sork on this system has been hampered by the limited amount of experimental data 
that is available. We discuss below one of the few theoretical attempts to understand 
the Bai_a;Ka;Bi03 system. The objective of this study is limited to understanding the 
insulating behaviour in BaBiOa system and its disappearence on doping with potassium. 



Rice and Wang interpret the BaBiOs system as a Pierls insulator which undergoes 
a semiconductor-metal transition on doping because of the loss of the commensurate 
Peierls condensation energy. They assume that the ground state at x = is a bond- 
ordered, CDW rather than an onsite CDW i.e., a Peierls insulator in which there is 



an alternation of bond lengths. Experiments indicate the possibility of p6| , |g, |2^, ^ 
alternating Bi-0 bond length in BaBiOa (though this bond-length alternation can equally 
well arise from a site diagonal CDW, as in the theories of Rice and Sneddon and Jurczek 
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and Rice). While considering the effect of doping, and the consequent disappearence of the 
insulating, bond ordered phase, they consider only the Bai_2,.Ki.Bi03, as the lead-doped 
system is diagonally disordered and a consideration of disorder is absolutely essential for 
its description. On the other hand, in the Bai-i^K^^BiOs, potassium replaces barium and 
hence the effect of disorder is probably less important. 

They start with the tight-binding model ||^, Sec. 1.3) involving the Bi6s conduc- 
tion band with (1 — x) electrons per unit cell in the Bai_^.K^.Bi03 system, hybridizing 
strongly with the three oxygen 2po- bands (with six electrons per unit cell ). The Hamil- 
tonian they use is: 

H='f Eiblb,, - ^ E - E VjM^Uab,. + H.c.) + C/2 u'^u,r (^^-18) 

j(T S jaS jS 

Here 6]^ creates an electron in a Bi6s orbital with energy 6^/2 (the zero of energy is 
set midway between the energies of Bi 6s and O 2p levels) at the site j of the reference 
cubic lattice (noncubic distortions are ignored); a]^ is a similar operator for the 2po- 
orbital (at energy -ep/2) and j + 6 denotes one of the six nearest neighbour sites of j; C 
denotes the spring constant for the stretching of the (j, j + 6) band, described by Ujj+s 
the deviation from the equilibrium separation 

Vjj+s = sign{6)(y - icUjj+s) (11.19) 

is the hopping matrix element between Bi and O where -7c is the derivative of Vjj^s 
with respect to u. The sign on the right hand side of Eq.(II.19) reflects the symmetry 
of the a bonds between Bi6s and 2p. Here the phonon couples to the hopping of the 
electron, whereas in the model introduced by Rice and Sneddon and Jurczek and Rice 
(Sec. II. 2. 2) the phonon couples to the charge difference on neighbouring Bi sites. In the 
case Ujj+s = for all j, the alloy is a metal and H is easily diagonalized to give two 
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cr bands and two non bonding dispersionless oxygen bands. In BaBiOs (i.e., at x = 0) 
the uppermost tight binding band (which is half filled) is unstable with respect to the 
commensurate Peierls distortion 

= (-l)'A/7e, foralH. (77.20) 

Here A is the amplitude of the dimerization and has the dimensions of energy. This 
dimerization leads to a splitting of both the original a bands and we get four subbands 
with energies £^A;,n=i- -4, given by 

E,,„ = ±[|' + 6(^2 + ^ [(i2V-A)2 + {y^ - t^^fellVW 

where = zV ^^^'^ ^ being the nearest-neighbour lattice vector. 

The Peierls gap is the gap between the upper two bands (or, equivalently, the gap 
between the two lower bands) and 

= + ^{y + A)2] V' - t-^ + 6(1^ - A)2] V'. (//■21) 

At half filling, the three lower subbands are filled and the topmost subband (from the 
splitting of the a* band) is empty. The total energy is then given by 

E^ot = 37VC(A/7e)2 + 2^^fe,i. (77.22) 

k 

The Peierls distortion parameter A is determined from the self-consistency condition 
obtained by setting dEtot/dls. = 0. A second derivative of Etot with respect to A gives the 
restoring force for small displacements of from its equilibrium value (as obtained 

from the self-consistency equation). 

According to this theory, then, the insulating gap of BaBiOa has to be equal to this 
Peierls gap. So the Peierls gap is fitted to the optical gap (2eV) of BaBiOa to obtain 
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estimates for the parameters of this theory. Using ep/2 = —1.1 eV and V = 2.2 eV (values 
taken from the band-structure results of Mattheiss and Hamann as discussed earlier) 
these authors estimate the Bi-0 spring constant C = 15.7 eV/A^ and the dimensionless 
electron-phonon coupling constant Agp = 2'y'^3'7iCV = 0.16. With these choices they 
get the bond-alteration amplitude at x = to be u ~ O.OSAwhich is quite close to the 
experimental value ||, jS^, |3T[ in BaBiOa. 

Next, an estimate can be made for the critical concentration xq of potassium at which 
the Peierls distortion vanishes and hence a semiconductor-metal transition takes place via 
the vanishing of A. This is easily done by partially emptying the second top-most (filled 
at x=0) subband and self-consistently calculating the value of A. xq turns out to be 
11%, considerably smaller than the experimental value of x = 0.40 for the Bai-^^K^^BiOs 
system. 

II. 2. 5 Local-Pair Superconductivity and Barium Bismuthates 

The local-pair model for superconductivity has been studied in detail by Micnas et. al 
||80| and several other people, mainly in the context of the possiblity of bipolaron super- 
conductivity |]85| , |90[| for a variety of systems (like Ti407) which show a CDW ordering. 
They studied the negative-f/, extended-Hubbard model in the largef/ limit. We do not 
review here the extensive literature on negative-?/ models which are not directly related 
to the bismuthates. 

The possible application of this local-pairing idea to the barium bismuthates has been 
pointed out by Varma [|18[. He proposed that the barium bismuthates be modelled by a 
negative-f/, extended-Hubbard model with strong onsite interaction U. In this limit, the 
model becomes identical to the local-pair model as mentioned above. He then noted that 
the phase diagram of the negative-f/, extended-Hubbard model is qualitatively similar to 
that of the barium bismuthates. 
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In the early local-pairing papers the negative-?/ was attributed to the electron-phonon 
interaction (as in Rice and Sneddon); however, the local pairing and bipolaron mechanism 
requires an extraordinarily strong coupling between electrons and phonons to produce an 
adequately large and negative-f/. Subsequently Varma has pointed out the possiblity of 
an electronic origin for the negative-f/ and, in a recent article on local pair supercon- 
ductivity, this alternative has also been alluded to. 

The Hamiltonian studied is the extended Hubbard Hamiltonian: 

H = t J2 4aCja-^J2^i<^^i~'^ + ^ H ^i^^-jV- (11.23) 

<ij>(7 ia <ij>aa' 

In the local-pairing model one takes U to be larger than the band width, so the 
negative-f/, extended-Hubbard model (II.2) can be mapped onto the pseudospin model 
(II. 5) and the phase diagram obtained via mean-field theory. Typical phase diagrams are 
qualitatively similar to those observed for the lead- and potassium-doped bismuthates. 
There are three distinct phases: (pseudospin) antiferromagnetic (CDW for model II. 2), 
spin flopped (superconductor for (II. 2)), and nonordered (metallic for model II. 2). A 
first-order boundary separates the superconducting phase from the CDW phase; all other 
transitions are continuous. In the T — x plane the first-order boundary should appear as a 
region of two-phase coexistence (although reported earlier such a region has not been 
confirmed in later experiments with high-quality samples). The superconducting phase 
of this model consists of strongly bound, onsite-paired electrons. A condensation of these 
bosons gives rise to supeconductivity. The system is, of course, diamagnetic as the spin 
degree of freedom is suppressed by the attractive interaction. We will discuss this model 
later in Section III and show that superconductivity and a host of other phenomena in 
the barium bismuthates cannot be described in terms of local pairing. 

II.3 Other Related Work 
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There have been some other recent attempts to understand the superconducting phase 
and high superconducting Tc in the barium bismuthates. An attractive interaction be- 
tween doped carriers in a CDW background can be generated, in certain circumstances, 



by the exchange of CDW amphtude modes The local depression of the CDW am- 
plitude because of doping leads to an attractive interaction for the doped carriers and 
mediates pairing in fc-space. The superconducting transition is of the BCS type, with the 
phonon cutoff frequency ujo replaced by a much higher scale of energy (the frequency of 
the CDW amplitude mode). The formalism is analogous to its positive-t/ counterpart. 



namely, the spin-bag mechanism ||124|| for high-Tc Cuprates. 

Some recent work has tried to understand the semiconducting behaviour of the 
Bai^i^K^^BiOs system starting from the model of Jurczek and Rice (Jurczek and Rice ||79| , 
see Sec. II. 2. 2). Their motivation was to estimate the various parameters of this model 
from the known experimental results on Bai_a;Ka;Bi03. 

They use the Jurczek- Rice Hamiltonian (II. 3) and assume that this system is in a CDW 
ordered state in its semiconducting phase. There are two types of sites in the system: Bi'^+ 
and Bi^+ at half filling with concentrations C3+ = C5+ = 0.5 each. Doping by potassium 
changes them to (1 — x)/2 and (1 + x)/2, respectively. In this picture the origin of the 
semiconducting phase is quite simple: The Bi^+ and Bi^+ sites give rise to the lower and 
upper bands containing 2c3+ and 2c5+ states and if there is a gap between them, then the 
Fermi level is in the gap as the filling of the lower band is 2c3+ = 1 — x. The system is then 
an insulator. Manh,, et al. calculate the band structure of the system described above 
where they account for the disorder due to potassium doping in the usual Bethe-Peierls 
approximation and the interaction is treated in the mean-field approximation. 

The local density of states is calculated by these authors in the alloy approximation: 
They use a continued fraction representation for one-particle Green function with three 
unknown parameters in the theory: onsite attraction U, nearest-neighbour Coulomb re- 
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pulsion V, and the effective electron-phonon coupling constant Aep (= gNlEp), where g 
is the electron-phonon coupling strength (Eq. (II. 3)) and N{Ef) is the DOS at the Fermi 
level. All the energies are scaled by the band width (the value used by the authors is 
the width of the a-antibonding band in Bai-xK^ BiOa, which is nearly 6 eV from band 
structure calculations). These authors try to fit the experimental values of the optical gap 
Eg at X = where Eg = 2eV and at x = 0.4 where, Eg = 0, (i.e., where Bai-^^K^^BiOs 
becomes metallic). With the gap in the DOS obtained by them at those values of x, they 
conclude that for —3eV < U < —2eV, and Aep < 0.7 produce a good fit for the optical 
gap. 

All the theoretical studies discussed so far are based on models where the onsite. 
Coulomb interaction is either suppressed (Rice-Sneddon case), or the onsite Ueff is either 
very small or even negative and attractive, and the nearest-neighbour Coulomb repulsion 



produces a CDW state. There is, however, a suggestion |^ that these systems are 
positive-?/ systems and their insulating phase is be a spin-Peierls insulator, melting at 
a high temperature into a resonating-valence bond liquid: One starts with a strongly 
coupled spin and lattice system plus a Heisenberg-Hubbard model for the spin and charge 
degrees of freedom. Now such a spin-Peierls state involving spins at each 31^+ ion will 
lead to a modulation in the Bi-Bi distance. As far as we know there is no indication of 
an alternating Bi-Bi distance in the structural work on these systems. We do not present 
a detailed discussion of this model, in this review. 

There is some work by Jin et al. [0] where a molecular dynamics simulation of phonon 
density of states is used to study the superconductivity in barium postassium bismuthate 
within a strong coupling Eliashberg formalism. These authors used a moderate electron- 
phonon coupling strength (A ~ 1) and assumed that strongth electron-phonon coupling 
exists even at high phonon energies (30-60meV). The Tcthat they could reach are fairly 
high. 
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II. 4 Discussion 



Most of the theoretical work on the bismuthates starts with: the negative-f/, extended 
Hubbard model to account for the experimental observation P, ^ [19| of a CDW phase 
close to half filling and the occurrence of superconductivity on doping. The diamagnetic 
behaviour of these systems can also be explained by using this model since the attractive 
interaction suppresses the free-electron moments. 

The theory of Sneddon and Rice was aimed at motivating the attractive U model 
for the barium bismuthate system. They did not try to obtain experimental fits from 
this model: A qualitative description of the pairing was their objective. In the model of 
Jurczek and Rice the more realistic case of finite band-width and disorder was discussed. 
The important result obtained from this work is that the CDW can survive rather large 
doping. The Tc obtained from this model was low for realistic electronic DOS at Ep- 
Their model produces a pseudo-gap at the Fermi surface, instead of a real gap as seen 
in the experiments, even in the insulating side and does not account for the observed 
transport properties in the insulating regime. 

The model of Yoshioka and Fukuyama focusses attention on the insulator-superconduc- 
tor transition; The connection between their "alloy" and the real BaPbi-jjEia-Oa system 
is rather far-fetched. The phase diagram obtained by them is symmetrical about x and 
1 — X, which is not the case for the real system. The theory of metal insulator transition 
in potassium doped system put forth by Rice and Wang [^, although interesting, is not 
quite satisfactory as it does not consider the effect of Coulomb interactions. Also the 
eventual development of superconductivity at low temperatures in the metallic regime of 
these systems is completely ignored. They obtaine a semiconductor-metal transition at a 
very low value of K concentration (11%) compared to the experimental value (~ 40%). 

The local pairing mechanism for the superconductivity in these systems essentially 



52 



involving the largef/ limit of the negative U, extended Hubbard model, can not produce 
reliable experimental fits of both the superconducting and the normal properties. It 
produces too small a superconducting coherence length and too large an optical gap. The 
theory of charge-bag can produce a superconducting Tc that is large but does not account 
for the observed unusual transport properties in the CDW or semiconducting phases of 
these systems. The positive U model has not so far been used to calculate the normal 
and superconducting properties. 

The unusual experimental observation of the wide separation between the optical 
and transport gaps has not been addressed by any of the theories we discussed. We 
find that this two-gap feature, coupled with other experimental observations, both in the 
superconducting and the normal phase of the bismuthates, puts quite stringent constraints 
on the models for these systems. We have recently explained this experimental observation 
and developed a theory for the bismuthate systems [ill which can explain many of 
the experimental observations in these systems on the basis of the negative-f/, extended- 
Hubbard model; we discuss this in the following sections. 

The negative-f/ generated by an electron-phonon coupling is retarded (and hence cut 
off in energy), whereas the negative-t/ generated by electronic processes acts over the 
entire non-interacting band- width. This makes the two mechanisms quite different, with 
different experimental manifestations. This is another major issue we focus on We argue 
that the two well-separated gaps can be understood only if the attraction is electronic in 
origin and suggest experiments to verify this assertions. We also calculate various normal- 
state and superconducting properties from the negative-t/ extended Hubbard model and 
the values of the parameters of this model that best fit experimental data. 

Though there are strong phenomenological grounds (Sec. II. 1) for using the negative-?/ 
extended Hubbard model for the bismuthates, there is no clear derivation of the negative- 
U. In our work, it is used as a parameter of the theory. In Sec. HI we discuss the large-t/ 
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limit of this model where a mapping onto a pseudospin, antiferromagnetic, anisotropic 
Heisenberg model (II. 5) allows us to investigate the spin model. We discuss briefly the 
mean-field theory for this model and obtain various correlation functions in the random 
phase approximation (RPA). We also suggest a novel mechanism of charge transport via 
pseudo-spin exchange. We discuss the limitations of the large— f/ limit and motivate the 
more interesting and appropriate limit of intermediate U which we take up in Section IV. 



III. The Negative-?/, Extended-Hubbard Model: Large-?/ Limit 
III.l: Summary of Results 

In Sec. II we motivated the use of the negative-?/, extended-Hubbard model (II. 2) for 
the barium bismuthates. It has been suggested [ITBI that, for a preliminary understanding. 



these bismuthates be modelled by the simpler Hamiltonian (II. 5), since theoretical studies 
of this model are far easier than those of model (II. 2). Model (II. 5), an anisotropic, pseu- 
dospin S = 1/2, Heisenberg antiferromagnet in a pseudomagnetic field, is the large— f/ 
limit {U ^ zt) of the negative-f/ (i.e., onsite attraction) Hubbard model (II.2) with 
nearest-neighbour repulsion V. (We follow the convention that U and V are always posi- 
tive; the sign of the interaction is displayed when needed.) It has been studied via differ- 
ent approximations by many workers p, |79|, 0] . Here we examine the viability of 



this limit for the barium bismuthates by doing detailed calculations for the phase diagram 
of this model and the properties - thermodynamic functions, correlation functions, and 
response functions - of its metallic, semiconducting, and superconducting phases 
and comparing them with experimental observations. Our overall conclusions are: (a) On 
the positive side, even the anisotropic-Heisenberg limit of the extended-Hubbard model 
leads to various results which are either in qualitative accord with experiments or which 
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suggest new experiments for the elucidation of the properties of these alloys, (b) However, 
it is not good enough to obtain many important properties of these systems; its most seri- 
ous limitation is that it leads to coherence lengths in the superconducting phase that are 
far too small compared to those obtained experimentally (~ 70 — 80 A). We are thus led 
to use the intermediate— [/ limit of model (II. 2), which we use in Sec. IV to obtain better 
agreement with experiments. We summarize our large-?/ results below: 

(1) In the simplest, mean-field approximation, the anisotropic, S = 1/2, Heisenberg anti- 

ferromagnet exhibits, typically, three phases, which, in the language of the extended- 
Hubbard model, are (seeSec.III for the correspondences between the two models): a 
nonordered metal (NO), a charge-ordered semiconductor (with a CDW), and an s- 
wave superconductor (SC). Phase transitions between all excepting superconducting 
and CDW phases are continuous (second order); the superconductor-CDW transi- 
tion is first order. Slight perturbations of the model (say by the inclusion of next- 
nearest-neighbour couplings) can yield an intermediate phase that exhibits both 
superconducting and CDW ordering |9^. At special values of coupling constants, 
bicritical, tetracritical, and tricritical points or critical end points can be obtained. 

(2) A random-phase-approximation (RPA) calculation of the excitation spectra in these 

phases shows the existence of a gap (~ V) in their semiconducting, CDW phase, 
which is different from the gap (~ U) for the breaking of pairs in the superconducting 
phase, which should appear in the optical spectrum. 

(3) It is useful to think of the phases in the following way: Bi^+ (6s^ configuration) 

corresponds to a real-space pair of two electrons and constitutes a charge-2e, hard- 
core boson. The CDW phase is a crystal of these bosons, the superconducting phase 
a superfluid resulting from their bose condensation, and the nonordered-metallic 
phase their normal-liquid phase. 
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Charge conduction occurs by the motion of these charged bosons, or equivalently, by 
the exchange of pairs of Bi^+ and Bi^"*" states. For this conduction mechanism the 
current operator has been constructed and thence, using the RPA, the frequency- 
dependent conductivity <j{uj) can be obtained [|10]- It is shown that: 

In the nonordered-metallic phase the conductivity ct(ci;) has a Drude form, i.e., cr(u;) = 
n(2e)^r/m*(l — iur), with r a phenomenological transport hfetime that accounts 
for scattering, m* the effective mass, and n the number of paired carriers; in the 
grand-canonical ensemble, i.e., with fixed chemical potential, goes to zero exponen- 
tially as the temperature goes to zero; with fixed (i.e., at fixed filling), this means 
that the metallic phase does not extend down to zero temperature at any filling (See 
Sec.IV). 

In the CDW semiconducting phase the conductivity shows an activated behaviour, 
i.e.,a{uj ^ 0) ^ exp {— A/ KbT), where A is the gap (measured from the chemical 
potential) in the spectrum of the hard-core bosons ( or, equivalently, the pseudospin- 
wave excitation spectrum ) in the CDW phase; this is in accord with experiments. 

Exactly as in the case of the bismuthates, the CDW phase of the model exhibits 
two gaps - a charge gap associated with the breaking of Bi'^+ pairs, which costs an 
energy~ U and another ( A ~ V^, the nearest-neighbour repulsion) associated with 
the activated behaviour of (t(c<j — > 0). We calculated the optical conductivity in 
the single- and two-particle channels using the Kubo formula and the RPA in the 
large-?/ limit. The result is shown in Fig. III. 4. It is interesting to note that even at 



a simple level of approximation ||T0|, the peak in the two-particle channel (at 0.24 



eV) is down by a factor of 9 compared to that in the one-particle channel (at 2 eV). 

The metallic phase of this model is not a conventional Fermi liquid, but rather a 
charged, hard-core, bose liquid, and must have exotic properties. For example, the 
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specific heat is certainly not a linear function of the temperature T. (A crossover to 
a conventional Fermi liquid occurs as U decreases.) 



We now show explicitly the calculations that lead to the above results. The calculations 
are described in two parts: The first part (III.2.1) is devoted mainly to to the mean-field 
theory of the anisotropic, pseudospin Heisenberg model (II. 5). In the following section 
(III. 2. 2) we show the calculations that lead to the excitation spectra in RPA in each 
of these mean field phases. We use the RPA excitation spectra then to obtain various 
correlation functions. We discuss the novel conductivity mechanism that we suggest for 
the bismuthates in their insulating phase and calculate the corresponding current operator 
in the large— f/ limit. 

III. 2 Calculations: An Outline 

III. 2.1 Mean-Field Phase Diagram 

The phase diagram of the pseudospin model (II. 5) is well known at the level of the 



so we 



mean-field approximation; however, some recent treatments are misleading [|80|, | 
outline briefly how these phase diagrams are obtained 0. 

The order parameters chosen are < S'f > and < > where i is the site index and 
can have only two labels A and B (indicating the two sublattices ). Also the rotational 
symmetry in pseudospin space about the 2;-axis (the magnetic field picks out the z direc- 
tion) allows one to choose the spin to be in the xz plane, i.e., < Sf >= 0, without any 
loss of generality. A straightforward mean-field decoupling scheme is used to reduce the 
Hamiltonian (II. 5) to the following site-diagonal mean-field Hamiltonian: 



nmf = J}2s!<s^> + }: st{K <s;> -b). 

<ij> <ij> 
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This Hamiltonian is diagonalized by the following eigenfunctions 



\z+ >= COS — I T> + sm — I l> 

and 

\i- >= COS— I i> - sm— I T>, 

where 

tan^, = Jo < > /{B - < S] >) 

and I |> and | |> are the eigenvetors of 5*^ operator with eigenvalues +1/2 and -1/2, re- 
spectively. Here Jo = and Kq = zK and [3 = (KbT)'^. The corresponding eigenvalues 
are 

Af = -Bt Rj 

with 

= (Jo < >y + {B-Ko< S; >y. 
The mean-field free energy per site is 

/ = -^[log(2cosh(/3i?^)) +log(2cosh(/3i?B))] 
+Jo < SI >< S% > ~Ko < S'^ >< S'b > -B. 

(This can also be obtained by using Peierls or Gibbs-Bogoliubov variational method. 
P7| Finally the self-consistency equations follow easily by minimising the free energy with 
respect to the order parameters: 

< 5f >= (Jo < 5J > /2Rj) tanh(/?i?j). (///.I) 

< >= [{B -Ko< SI >)/2Rj] tanh(/?i?^); (^^^-2) 
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The subscripts i and j label nearest-neighbour sites; i.e., they belong to the two 
different sublattices A and B (or B and A), respectively. If more than one solution is 
obtained for these self-consistency equations, then the one which yields the lowest value 
of the free energy is picked. 

Note that we can work with either fixed pseudomagnetic field (the analogue of the 
grand-canonical or fixed-chemical-potential ensemble for the original, extended-Hubbard 
model) or fixed < 5*^ > (the analogue of the canonical or fixed- filling ensemble); the 
results obtained from both approaches have to be the same by virtue of the equivalence 
of the two ensembles (see Note 2) . 

The phase diagram in both canonical and grand-canonical ensembles have been ob- 
tained. A typical phase diagram that results from this mean-field calculation is shown, 
in Fig. III. la, in the B — T plane and, in Fig.III.lb, in the < 5*^ > — T plane. (Pa- 
rameters for this phase diagram are: Jq = 1, -f^o — 1-301.) There are three phases: 
a nonordered metal at large B or T (the pseudospin-model analogue is a paramag- 
net with some alignment of the pseudospins by virtue of the finite magnetic field), a 
CDW semiconductor (the pseudospin-model analogue: a ferrimagnet (antiferromagnet at 
B — 0)), and an s-wave superconductor (pseudospin-model analog: a spin- flopped phase, 
with nonzero x component of pseudospins). The pseudospin-model analogues can be 
obtained easily by consulting Table II. 2. In terms of the hard-core bosons introduced 
above, these three phases correspond, respectively, to a bose liquid, a bose crystal, and a 
bose-condensed superfiuid (note 3) phase. 

Full hues in Fig. III. la indicate continuous, second-order transitions, whereas the 
dashed line in Fig. III. la indicates a first-order phase boundary; this boundary shows up 
as a region of two-phase coexistence if < 5"^ > (i.e., the filling) is held fixed (Fig.III.lb). 
If the temperature and filling are such that the system lies in this two-phase region, then, 
in equilibrium, the model system phase separates into two regions - one superconduct- 
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ing and the other semiconducting - with one interface between them. This two-phase 
region is not a new mixed phase with both superconducting and CDW order. To obtain 
such a phase, next-nearest-neighbour interactions are required in the pseudospin model 
(II. 5) (see below). Note that at no value of the filling does the metallic phase extend 
down to T = 0. (This manifests itself in an exponential depletion of the carrier den- 
sity if one tries to approach T = at fixed B .) The two second-order lines meet the 
first-order boundary (Fig. III. la) at a bicritical point. The positions of these second-order 



lines can be obtained analytically as described in Appendix B. As Varma has noted p!8|] 
in this pseudospin-model limit the maximum Tc for the superconducting phase occurs at 
the value of filling for which the metal-semiconductor transition first appears, i.e., at the 
filling at which the bicritical point occurs. (Strictly speaking, this maximum occurs at a 
slightly higher temperature. The two critical lines come in tangentially at the bicritical 
point, but this becomes evident only from a renormalization-group calculation.) 



It is well known from the work of Liu and Fisher and Matsuda and Tsuneto 
that slight modifications of the pseudospin model (II. 5) can lead to more complicated 
phase diagrams, at least at the level of mean-field theory. In particular, if we allow for 
next-nearest-neighbour interactions of the form JijSi-Sj, where << ij » in- 

dicates that i and j are next-nearest neighbours, then, depending on the relative sizes 
of the different interactions in the model, phase diagrams of the sort shown schemati- 
cally in Fig.III.2 can be obtained. (Such next-nearest-neighbour interactions could arise 
from next-nearest-neighbour repulsion V or hopping t' in the original extended-Hubbard 
model.) With such couplings it is possible to get a bona fide intermediate phase in which 
both CDW and superconducting order are present (Fig.III.2). As in Fig. III. 1 solid lines 
denote continuous transitions and dashed ones first-order transitions. Note that a variety 
of multicritical points can be obtained now [|98|, |9^: bicritical, tricritical, and tetracriti- 



cal points, and also a critical end point. We mention these multicritical points here and 
discuss them briefly in Appendix D mainly for the sake of completeness. It is quite likely 

60 



of course that the precise multicritical behavior that obtains in the barium bismuthates 
depends sensitively on effects that are not included in our model, such as disorder, long- 
range Coulomb interactions and deviations from a cubic crystal structure, etc. Also, it will 
not be very easy to determine the nature of such multicritical behavior experimentally, 
for that will require very accurate control of both the temperature and the concentration 
of the dopant. 

The phase boundary between the CDW and nonordered metallic phases belongs to 
the three-dimensional Ising universality class and the superconductor-nonordered metal 
transition to the three-dimensional XY class. The transitions from the intermediate 
phase (Fig. III. 2) to the CDW and superconductor phases are, respectively, in the three- 
dimensional XY and Ising universality classes. In the pseudospin model, the coherence 
length in the superconducting phase is very small (see Appendix F), so apparent mean- 
field exponents will not be seen here as they might well be in the barium bismuthates, in 
which coherence lengths are of the order of 80 A. 

III. 2. 2 Excitation spectra, Correlation function and the Optical Conductivity 
in the Random-Phase Approximation 

The excitation spectra in all the phases of the pseudospin model can be obtained 
in a straightforward manner by using the equation-of-motion method followed by the 



random- phase approximation ||10(]|| as shown in Appendix C. 



We use the Hubbard operators ||100|| P^^ = \ia >< i(3\, where the state \ia > refers to 
the states |+ > or |— > at the site i, in the mean- field Hamiltonian described above. The 
double-time retarded Green functions for these operators are denoted by G^£^g(t — t') or, 
equivalently, << P^p{t)\PJ^g{t') ». The Green function is defined in the standard way 
as 
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where [P,Q]_ denotes a commutator between the operators P and Q. 

The equations of motion for the Green functions are easily obtained in the random- 



phase approximation The poles of these Green functions can be obtained simply 

from the solutions to the equations of motions. From these poles we can obtain the 
excitation spectra in the various phases of Fig.III.l. The details are discussed in the 
Appendix C. 

In Fig. III. 3 we show the pseudospin-wave excitation spectra in the three phases 
(respectively, CDW, superconductor, and nonordered metal) of the phase diagram of 
Fig. III. la. The spectra are drawn along the symmetry directions of a cubic Brillouin- 
zone. Data for these spectra are given in the figure captions. Note that there is a gap as 
expected in the CDW (i.e., ferrimagnetic in the pseudospin language) phase. This gap 
is obtained from the dispersion relation in the CDW phase: 



uj^ = 2B- Kom ± [K^ + (m^ - r]^)J\q)]. {1 1 1 A) 

Here 77 =< S'f > — < > is the staggered magnetisation. The ± sign indicates that 
there are two branches in the spectrum, as should be the case in an antiferromagnetically 
ordered spin system. There is a gap in the spectrum of the non-ordered metallic phase also 
at any finite magnetic field which is merely the chemical-potential gap of the charged bose 
liquid. (In the pseudospin language this is the gap due to the aligning pseudomagnetic 
field.) From this expression it is easy to see that the gap at T = in the CDW phase is 

A = [K^ - J'(0)]^/^ (111.5) 

Where 
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^(q) = ^07q, 



7q = -l]exp(iq.a), 

^ a 

with {a} is any one of the nearest-neighbour lattice vectors, 
and 

m = n-l^<S-> + <Sj> 

The gap A vanishes ( as does the CDW phase) as the anisotropy goes to zero, i.e., as K 
starts to become equal to J. This is the transport gap Ea that shows up in measurements 
of the dc conductivity in the CDW phase see below. 

The excitation spectrum for the metallic phase also shows a gap. However, this does 
not lead to activated behaviour in the conductivity: This gap is merely the chemical- 
potential gap of the charged bose liquid that constitutes the metallic phase and is a 
manifestation of the finite charge susceptibility a,t q — and T > (see Note 1). 

Various correlation functions can be obtained from the Green functions that have been 
calculated. For example, in the nonordered phase 

< S+{q,uj)S-{q,uj) >= ^{n - l)/{uj - e(q)), 

where 

e(q) = 2S - (n - 1)7^0 - (n - 1) J(q) = P'' log(n/(2 - n)) - (n - 1) J(q). 
Similarly in the CDW phase 



1 [u ~ 2B - 2Ko{m + ri){m + 7])] 



(777.6) 
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It is easy to check that at = these correlation functions diverge as g ^ at 
the critical and multicritical points of the phase diagrams of Figs. III.l and III. 2. Here 
this divergence is characterized, of course, by mean-field exponents; corrections to these 
exponents and to the scaling forms for various thermodynamic and correlation functions 
can be obtained from renormalization-group calculations as described in Appendix D. 

Given the Green functions that have been calculated, we can easily obtain the 
frequency-dependent conductivity by noting that the only mode of conduction in the 
limit of large U is that in which the charged boson comprising Bi^+ moves or, equiva- 
lently, pairs of Bi'^+ and Bi^+ sites are exchanged. The current operator for this mode of 
conduction, (Appendix E and is given by: 

3 = -t{2e)Jj2iSt-Si-H.C.)B.. 

<ij> 

By using the Kubo formula and an RPA- decoupling scheme we can evaluate the con- 
ductivity in all the phases ( all effects of scattering are contained in a phenomenological 
scattering time r) . The result of such a calculation for the metallic phase is of the Drude 
form 

, , n(2e)^T/m* 

(1 — lUT) 

where n is the number of charge-2e carriers and m* is the effective mass. 

For the CDW phase, as — > we get a{u! — > 0) ~ exp(— A/ff^T). (This is the anisotropy 
gap in the pseudospin-wave excitation spectrum in the antiferromagnetic phase.) 

III. 3 Remarks on the viability of the Large— Umit for the bis- 
mut hates 

The phase diagram and the other results obtained above for the pseudospin model have 
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many features in common with the observed properties of the barium bismuthates (Sec. 
III. 2.1). However, there are serious hmitations which, in our opinion, make this large-?/ 
model not viable. The negative-f/ Hubbard model may be a more appropriate limit for 
the barium bismuthates and may offer a more realistic description of these systems if we 
choose intermediate values of U. 

The most serious limitation of the large-f/ pseudospin model is that the coherence 
length in its superconducting phase is about one lattice spacing. This is obvious since the 
model corresponds to hard-core, charge-2e bosons which are onsite-paired electrons. An 
explicit calculation of this coherence length is given in the Appendix F. In contrast, the 
experimentally measured coherence length is much larger (~ 80 A). 

One of the attractive aspects of the large-[/ pseudospin model is that it yields the 
two gaps seen experimentally (see Sees. I and II) - an optical gap ^ U and a transport 
gap ~ 2zV. From the experimental data for BaBiOa one would therefore have to infer 



that U ~ 2.0 eV and 2zV ~ 0.24 eV. However, from the DOS it can be estimated 
that (see Note 4) zt = O.SeV, so the large-[/ assumption (i.e., U ^ 2zt) is clearly not 
true. Also the charge contrast (relative to uniform Bi^"*") implicit in the CDW phase of 
the large-f/ model is too large compared to its value in these bismuthates. (The charge 
contrast has not been measured precisely in any experiment. However, the value of the 
charge contrast that can be inferred from the relative sizes of the enlarged and contracted 
oxygen octahedra is much smaller than ±1, which the large-U model predicts). 

Both the shortcomings discussed above show that, if the negative-f/ model is at all 
viable for these bismuthates, it is so only if the large-f/ limit is abandoned. Thus we must 
look at the intermediate- and small-f/ limits, which we discuss next in Sec. IV. 
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IV. The Negative-^/, Extended-Hubbard Model At Intermediate 

U 



IV. 1 Summary of Results 

In this Section we continue our discussion of the viabihty of modehing barium bis- 
muthates using the negative-?/, extended-Hubbard model (with a repulsive interaction V 
between nearest-neighbour sites) in the intermediate-U case. We obtain the phase dia- 
gram of this model via a mean-field theory and the compute the effect of fiuctuations at 
T = in the RPA. The details of our calculations are discussed elsewhere [1^, |TT|]. 
We discuss only the principal results here. Furthermore, we comment on the feasibilty of 
using the negative-t/, extended-Hubbard model for these bismuthates by fitting a variety 
of experimentally measured quantities, such as the optical and transport-activation gaps 
in the semiconducting phase, the superconducting Tc and coherence length j^, |10|, |lT . 

Our principal results for the intermediate-?/ are: 

(1) A phase diagram (Figs. IV. la and IV.lb) for the intermediate- f/ case, extended- 
Hubbard model by using a mean-field approximation which is can distinguish be- 
tween a variety of phases: a semiconducting, charge- density- wave (CDW) phase, 
superconducting phases with s-wave, extended-s-wave, or both of these types of 
ordering, intermediate phases with both superconducting and CDW order, and a 
nonordered, metallic phase. This phase diagram has the same topology as that of 
the \axge-U limit (Sec. HI). Furthermore, it is in qualitative accord with experimen- 
tal phase diagrams in so far as it exhibits three distinct thermodynamically stable 
phases: a CDW phase, a singlet superconducting phase, and a nonordered metallic 
phase. Phase transitions between all excepting CDW and superconducting phases 
are continuous; the superconductor-CDW transition is first order. 
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If the negative-f/ is electronic in origin (see point (6) below), then in the CDW phase 
there are two-particle and two-hole bound states, i.e., Cooperons, with energies in 
the gap of the two-particle spectrum, (Fig.IV.2b). We have calculated the pairing 



susceptibility within the RPA in the particle-particle channel |TOl to obtain these 
bound-state (Cooperon) energies in the CDW phase. In the limit of large U, this 
energy evolves smoothly into the anisotropy gap in the pseudospin-wave spectrum 
(Sec. III). The bound-state energy moves continuously towards the chemical potential 
as V approaches zero. 

Charge transport in the CDW semiconducting phase is dominated by the motion 
of these Cooperons. The large-f/ limit of this process reproduces the conductivity 
mechanism in the pseudospin model (Sec. III). Within the RPA the conductivity 
shows an activated behaviour with an activation energy equal to the energy of 
the particle-particle bound state. This is the transport gap. It is thus clearly 
different from the optical gap. For the latter is equal to the gap in the single-particle 
spectrum, i.e., it is twice the CDW gap parameter in the small-t/ regime. At large 
values of U, this gap is ~ f/; note that this is not the gap in the pseudospin-wave 
spectrum, but is the energy required to break a pair. 

At large U the pairs are tightly bound in real space and can be thought of as a 
quantum lattice gas of hard-core bosons; at small U k-space pairing is the proper 
description of these pairs. This calculation yields a natural interpolation between 
these two limits at T = 0. 

The conductivity in the nonordered, metallic phase is of the Drude form. At large U 
the carriers are charge-2e, hard-core bosons (Sec. Ill); at small U the metal is, of 
course, a conventional Fermi liquid. 

For the above explanation for the two gaps to be valid, the attractive, onsite inter- 
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action must be electronic in origin: An attraction mediated by phonons is retarded 
and operates only over some characteristic phonon frequency. Such an interaction 
cannot produce the necessary binding to form the two-particle bound state as the 
CDW gap is ~ lev, and any phonon-mediated attraction must be cut off well below 
this energy range. 

The remaining part of this Section is organized as follows: In Section IV. 2 we describe 
our mean-field and RPA calculations in the intermediate-f/ range of the negative-f/ Hub- 
bard model very briefly. In Sec.IV.3 we conclude with a discussion of the strengths and 
weaknesses of the model, suggest future directions of study, and propose some experiments 
which can be used to verify our theory. 

IV. 2 Outline of the Calculation 

IV. 2.1 Mean-Field Theory and The Phase Diagram 

The one-band, negative-?/, extended-Hubbard Hamiltonian, defined on a simple cubic 
lattice, is given by 



<ij><7 icr <ij>(ja' icr 

Though the qualitative features of the phase diagram of model (IV. 1) are known 
||101|| , some earlier mean- field treatments ||101|| of it contain some errors (see Note 5). The 



qualitative features of the phase diagram of model (IV. 1) remain the same as those of the 
large-t/ model; however, there are quantitative differences. 

The phase diagram (Figs. IV. la and IV. lb) is obtained by using a mean- field approx- 
imation which is capable of distinguishing between a variety of phases: a semiconducting, 
charge- density- wave (CDW) phase, superconducting phases with s-wave, extended-s-wave. 
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or both of these types of ordering, intermediate phases with both superconducting and 
CDW order, and a nonordered, metaUic phase. The order parameters used are: 

< >= 

< CiaCja >= Co, 

and 

< n,^ I + |e^^-"-% {IV.2) 

where Q = 7r(l, 1, 1), be is the CDW order parameter, bs and be are, respectively, onsite and 
extended-singlet-superconducting order parameters, and Cq a self-consistent contribution 
to the band width. Only translationally invariant order parameters are considered. A 
standard, Hartree-Fock decoupling, yields the mean-field Hamiltonian (we take the order 
parameters to be real) 



k<T k(T k 

+N[Ubl - Vzbl + -^^Ura + ^-\?eU^ + V zCl\ . {IV. 'i) 

Here N is the total number of lattice sites, Cko- is the Fourier transform of Cjo-, Um = 
U — 2zV, Up = U + 2zV, p, = /I + Umn, ek = zi^\^; 7k = ^ X^a^^Pl^^-^)) where {a} 
are the nearest neighbour lattice vectors; Ac = Upbe/2 and A5(k) = Ubg — zVbe^k are, 
respectively, the gap parameters in the CDW and superconducting phases. 

This mean-field Hamiltonian (IV. 3) can be easily diagonalized and the four eigenvalues 
are ± E^, where 
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and 

With these four eigenvalues the mean-field free energy per site becomes 



In 4 cosh + In 4 cosh 



1 

+Uhl - Vzbl + -n\m + — + zVCl + Ji{n - 1) 

4 Up 

The order parameters are determined from the self-consistency equations that follow 
from a minimization of this free energy. The prime over summation indicates a reduced 
brillouin zone. A numerical, iteration method is used to solve the self-consistency equa- 
tions for the order parameters for different values of U/zt and V/zt (z, the number of 
nearest-neighbour sites, is 6 in three dimensions) and at various values of filling 5 = n — \. 

The phase diagrams in both filling-temperature and chemical potential-temperature 
planes are shown in Figs. IV. la and IV. lb. The topology of the phase diagram remains 
the same as in the large-t/ case, but the numbers obtained (for the positions of phase 
boundaries, etc.) are quite different. There are three thermodynamically stable phases: 
A CDW phase characterized by Ac 7^ and = 6e = 0; separated from this CDW 
phase by a first-order line in the 11 — T plane is the superconducting phase (characterized 
by 6s 7^ 0, 6e 7^ and Ac = 0); and then there is a nonordered phase (all the order 
parameters are zero in this phase) at high temperature (at 5 = \ this phase goes down to 
T — Q). The phase transitions between the ordered phases and the nonordered phase are 
all continuous. The two lines of continuous transitions meet at a bicritical point. The first- 
order line in the grand-canonical (i.e., fixed /x. Fig. IV. lb) ensemble opens out into a region 
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of two-phase coexistence (CDW and superconducting phases) in the canonical (i.e., fixed 
6, Fig.IV.la) ensemble. Co is nonzero in all the phases. The hatched region in Fig.IV.la 
is the two-phase coexistence region. In this two-phase region (and in equilibrium) the 
system phase separates into semiconducting (CDW) and superconducting phases, with 
one interface between them. 

For V = and at half-filling (i.e., 5 = 0) in model (IV. 1) the CDW and supercon- 
ductor solutions of the self-consistency equations are connected by a pseudospin-rotation 
symmetry (see Note 6) of the Hamiltonian and hence are degenerate; their free ener- 
gies are equal and they coexist. Away from half filling this degeneracy is removed and 
the superconducting phase becomes the thermodynamically stable one: In the large-C/ 
limit this corresponds to the (pseudo-spin) antiferromagnetic phase giving way to the 
(pscudo) spin-flop phase in a (pseudo)magnetic field. The degeneracy between CDW and 
superconducting phases is also lifted when V becomes nonzero (the pseudospin analogue, 
described in Table II. 2, is that the pseudospin model of Sec.III becomes anisotropic for 
V 0): at half fiUing the CDW phase is the thermodynamically stable one; of course, 
sufficiently far away from half filling, the superconducting phase becomes stable. 

The gap equation for the CDW gap parameter is of the usual BCS form, from which 
the CDW-nonordered transition temperature follows. The gap equation is 



where 2W is the noninteracting band width. At half fiUing p, — and, if we take 
the density of states p{e) to be a constant, this becomes the usual BCS gap equa- 
tion. At T = 1/P — and with a constant DOS this gives the CDW gap parameter 
Ac^W/smh{2W/up). 

For the superconducting phase the calculations are also of the BCS type. The gap 
equation now involves the 2x2 matrix equation: 
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M 



1 ^ tanh(ek + fJ-) f U -Vz'j] 



iIV.5) 



2N^ (ek + /i) V -Vzji 

where 1 is the 2x2 unit matrix. The determinant of the matrix M should vanish at 
Tc. Note that, in this mean-field theory, the maximum superconducting occurs at the 
bicritical point. 

IV. 2. 2 The Cooperon Bound State In The Random Phase Approximation 

We showed in Sec. Ill that the experimental observation [0] of well-separated optical 
and transport gaps can be explained by identifying the optical gap with the pair-breaking 
energy and the transport gap with the pseudospin-wave gap in the CDW semiconduct- 
ing phase. This explanation remains viable even in case of the negative-U , extended- 
Hubbardmodel for intermediate V. 

As in the large-?/ case the CDW phase continues to show two well-separated gaps 
for intermediate values of U and V. The physics of these two gaps is as follows: Our 
mean-field theory yields quasiparticle valence and conduction bands with a gap separating 
them, (Fig.IV.2a). In k-space this gap opens up along the entire noninteracting Fermi 
surface. In the ground state, the valence band is full whereas the conduction band is 
empty (at half filling, i.e., for pure BaBiOs) and this system is insulating. Now consider 
putting two electrons into the conduction band. They still feel the residual attraction U 
(modified by appropriate coherence factors in general). Hence they form a bound state 
(a Cooperon, i.e., a charge-2e boson) in the gap in the two-particle excitation spectrum, 
(Fig.IV.2b). Two holes in the valence band form a similar bound state. This bound state 



forms provided the source of the attractive U is electronic [|1^] and not phononic since 
U must remain attractive over energy scales larger than the CDW gap Ac, which is 1 
eV, as measured from the Fermi level. Also a phononic mechanism cannot give rise to a 
binding of nearly 0.8 eV (required to fit experimental observations; see below). It is easy 
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to see that charge (and all other) transport in this CDW (semiconducting) phase will be 
dominated by these Cooperons. The thermal activation energy required for exciting these 
Cooperons is the gap for exciting two quasiparticles minus the binding energy of these 
Cooperons (i.e., the energy of the bound state measured from the chemical potential) 
(Fig.IV.2b). Hence it is lower than (half) the semiconducting (CDW) gap 2Ac. On 
the other hand, optical experiments, which excite single particles across the CDW gap, 
measure this larger gap, Ac. 

The existence of a two-particle bound state in the CDW phase is indicated by the 
appearence of a pole in the pair susceptibility in the Cooper channel (at wave vector 
q = 0) below the bottom of the continuum of the two-particle excitation spectrum (the 
two-particle continuum starts at 2 Ac (Fig.IV.2b). 

The (singlet) pair susceptibility in the Cooper channel is defined by 

/ + 00 ' 
dr exp(iz/r) ^ (T, [ck+qT,i' {r)c^H,j' (^)4'+gT,*(0)^-k'i,j (0)]) ' 
k,k' 

where represents time-ordering, i,j, etc. are band indices (in the CDW ground state 
there are two bands), the subscript p stands for pairing, and the expectation value is 
calculated with respect to the interacting ground state. We restrict ourselves to T=0. 

We show in Fig. IV. 3 the RPA ladder diagrams for the pair susceptibility Xp,ij;i'j'i^^ 
We have calculated the sum of these diagrams 0, |1T| in the presence of a CDW 
background and starting from a negative-f/, extended- Hubbard model at half filling. From 
this the two-particle (hole) excitation spectrum is obtained, and thence the location of 
the Cooperon bound state. It is shown that, as long as V (more precisely, A^), is small 
which is satisfied for the barium bismuthates because the A^ — leV), the bound state 
falls within the gap of the two-particle spectrum (Fig.IV.2b) and, for any nonzero V , 
this state lies above the chemical potential (all energies are measured from the chemical 
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potential). Hence, so long as V ^ 0, the CDW phase is not unstable with respect to the 
superconducting phase at half filling. 

We show that the resulting two-particle bound-state energy in the intermediate-[/ 
regime evolves into the pseudospin-wave gap of the large-U model. It is interesting to 
note that the RPA excitation spectrum calculated at intermediate values of U evolves 
into the RPA pseudospin-wave excitation calculated in the \aige-U limit (for it is not 
evident a priori that the two RPA's are equivalent in all respects). 

In the presence of a CDW background we have done the analogue of the Cooper- 
instability ||102|| calculation by following the t-matrix approach ||103|| . This requires the 
solution of a secular equation involving a 14x14 matrix ||^, |10|. Had it been possible to 
cast this equation into an equation of the type 1/A = $(z^) (where A is some effective, at- 
tractive interaction), then its solutions could have been portrayed as shown schematically 
in Fig.IV.4, i.e., as the intersections (empty circles) of the horizontal line {y = 1/X, where 
A is negative since the residual interaction is attractive) with the curves (thick lines) rep- 
resenting $(i^). The vertical lines (thin lines) are the unperturbed CDW eigen-energies 
starting from ±2Ac. We have taken the system to be finite system for the purposes of 
this figure, so the energy levels are discrete. As the system size goes to infinity, the eigen- 
energies of the interacting system which fall above 2Ac or below -2Ac form the two-particle 
and two-hole continua respectively, and correspond to scattering states; however, the two 
eigenenergies inside the gap remain separated from these continua, and correspond to the 
bound states (the one below the chemical potential is the two-hole bound state and the 
one above is the two-particle bound state). 

The energy of the bound state depends on U, V and zt. In order to pin down at least 
some of the parameter values, we proceed as follows: 

(1) Set Ac = leV, since the optical gap, which is twice the CDW gap parameter, is 2eV 
for BaBiOs. (We use the mean-field gap equation and a constant density of states which 
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makes /\c = zt/ smh{2zt / Up) , where Up = U + 2zV and zt is half the band width.) 
(2) The Cooperon bound-state energy (measured from the chemical potential) is set equal 
to the transport gap of BaBiOs, i.e., 0.24eV. The values of U (both screened and un- 
screened; see below) and V required to satisfy these two conditions for different choices of 
the values of zt are plotted in Figs. IV.5a and IV.5b. The main results obtained from the 
above calculations, relevant to the modelling of the barium bismuthates, are summarized 
below: 

(1) For the potassium-doped system the band width, close to the metal-semiconductor 



transition is roughly 1.6eV, as obtained from the DOS (Batlogg, et al |24|). Given 
this value we obtain, from Figs.IV.5a,b U = 1.9eV and 2zV = 0.13eV. Clearly then 
we must consider the intermediate-?/ and not the large-?/ case. 

(2) The two-particle bound-state energy goes to zero continuously as V approaches zero. 

At V=0 there is a gapless two-particle mode that signals the instability of the 
CDW state towards superconducting order. In Fig.IV.6 we show the values of the 
interaction-parameters obtained from the two conditions described above, for a fixed 
value of W{= zt), as a function of the bound-state energy. Note that the extended- 
Coulomb term V goes to zero as the energy of the bound state approaches the 
chemical potential. 

(3) The bound-state energy goes to 2zV as the band width approaches zero (Fig.IV.5b). 

Note that in the large-U limit the anisotropy gap is exactly equal to 2zV, a feature 
that is reproduced correctly here as ^ 0. This shows that the above results inter- 
polate smoothly between the small- and large-f/ results for the extended-Hubbard 
model at T = 0. 



We now compare our results with experiment als. However, we cannot naively use 
our estimates for the values of U and V to calculate the superconducting and CDW 
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transition temperatures for the following reasons: Our estimates of U and V have been 
made at the extreme insulating limit of BaBiOs, whereas superconductivity and the metal- 
semiconductor transition occur after significant doping ( 5 ~ 0.4 for the potassium-doped 
system). We must, therefore, try to incorporate the effects of metallic screening on the 
strength of the Coulomb interaction. In addition, doping has another important effect on 
these systems. The arguments used for an attractive, onsite interaction were based on 
nonlinear screening of the ionised configurations of the Bi atoms by the highly polarizable 
O octahedra. It is conceivable that, with increased doping, as the system moves towards 
the metallic side and the lattice distortions and the charge-disproportionated, highly 
ionized Bi configurations begin to disappear, the mechanism which leads to the attractive 
interaction gets suppressed rapidly. At this stage it is difficult to calculate these effects, 
as we do not have a first-principles theory for the negative value of U . 

In order to account for at least some of the physical processes that result from dop- 
ing, we have included the effect of metallic screening at the simplest level via the RPA. 
The corresponding RPA screening diagrams are shown in Figs. IV. 7a and IV. 7b (repeated 
excitations of particle-hole or particle-particle pairs, out of the Fermi sea ). Note that 
screening operates on both the particle-particle and particle-hole channels. Hence the ef- 
fective interactions responsible for CDW and superconductivity are both reduced. Using 
the standard RPA procedure to calculate the t- matrix we write the corresponding effective 
(screened) interaction as Vsc(q) = ^(q)/ [1 ~ ^(q)Xo(q)] where V is the bare interaction 
and Xo(q) is the appropriate bare susceptibility (particle-particle or particle-hole). Now 
the particle-hole susceptibility at wave vector Q is the same as particle-particle suscep- 
tibility at zero wave vector. Hence, in the simplest approximation, we can replace Xq by 
the density of states A'"(0) at the Fermi level and get the effective screened interaction 
governing the CDW or the superconducting transition. For A^(0) we use the experimental 
value at the metal-semiconductor boundary of the potassium-doped system. 
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By using the above procedure to account for screening, we plot the results for the 
superconducting and CDW transition temperatures in Figs. IV. 8a and IV. 8b calculted 
via the BCS formula discussed earlier (with both bare and screened interactions), as a 
function of the band width. The results for the phase boundaries obtained using this 
procedure are discussed in Appendix C and shown in Fig.IV.9. Note that we have used 
screening corrections only for the on-site Coulomb term. A calculation of the t-matrix for 
the momentum dependent V term is rather complicated and is not done here. Since the 
values of V we are interested in are quite low (Fig.IV.9) compared to U, it is reasonable 
to neglect the screening corrections to V. Indeed the dashed line in Fig.IV.9 shows the 
superconductor-nonordered phase boundary with V=0 (we show the results for two dif- 
ferent set of parameters values, and note that the parameter values that are appropriate 
for the bismuthates correspond to Fig.IV.9). This shows the negligible effect V has on Tc 
in the range of interactions that we are interested in. In passing we note that the effect 
of V on the superconducting Tc vanishes as 6 approaches zero |T0|. 



In any case, our results clearly show that the values obtained for the superconducting 
and CDW transition temperatures are too high compared to the experimental values. The 
effects of fluctuations neglected in our mean-field theory and the effects of disorder could 
reduce these transition temperatures by as much as a factor of 4, but probably not much 
more. It is likely that the reduction of the strength of attraction with doping (as described 
earlier) must be incorporated in our theory order to get better experimental fits for the 
Tc. (This especially for the superconducting transition temperature; the mean-field, CDW 
transition temperature is not orders of magnitude off from the experimental value, lending 
support to our hypothesis that, as doping increases, the negative-?/ is reduced.) 

In spite of this limitation of large transition temperatures, our theory yields quite a 
few interesting results. Firstly, the large difference between the transport gap the optical 
and gap can be understood from the existence of the two-particle bound state in the 
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CDW phase. The conductivity is given by cr(a;) = (ne*^r/m*)/(l — iuir), where m* is 
the effective mass of the charge carriers, i.e., the Cooperons and e* = 2e is their charge. 
n{T) exp(— E^^/XsT) is the number of charge-2e carriers at temperature T, where Ea 
is the transport gap. Hence (7(0; —> 0) shows an activated behaviour with an activation 
energy Ea- Secondly, this theory clearly reproduces some of the qualitative features of the 
phase diagram for the bismuthates (see the concluding discussion below). In addition, 
the value of 2A/KbTc is 3.5, and the system is diamagnetic at low temperatures in 
our model. We have taken the standard mean-field expression for the superconducting 
coherence length and used the observed values of the Fermi velocity to obtain a coherence 
length of nearly 30A. Since our calculation grossly overestimates the superconducting Tc, 
it is not surprising that our estimate for the coherence length is quite small compared to 
the experimental value. Nevertheless, it is much better than the extremely short coherence 
length predicted by the large-?/ (Sec. Ill) where the superconductivity arose from onsite 
paired electrons. Here, in the intermediate- U case, the pairing is in momentum space and 
the pairs are spatially extended. 
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IV. 3 Concluding Discussions 



In this and the preceding Sections we have investigated the negative-C/, extended- 
Hubbard model to model the barium bismuthate superconductors. In our study we have 
assumed that the mechanism that leads to the negative value of U is electronic rather than 
phonon mediated. At the end of studying both large- and intermediate-[/ ranges we have 
been successful in understanding some of these features associated with the bismuthate 
systems. 

Perhaps the most important aspect of our treatment is the detailed theory we have 
provided for the occurence of two gaps in the semiconducting phase of the barium bis- 
muthates. We have identified the transport gap as the energy of the two-particle bound 
state (Cooperon) and the optical gap as twice the CDW gap. This picture of the two 
gaps is consistent with the fact that the transport gap is not easily accessible to optical 
refiectance, photo-acoustic or photoconductivity studies, all of which involve single parti- 
cle excitations. Note that the above picture implies that the semiconducting phase of the 
bismuthates is unique, with transport properties dominated by charge-2e bosons. 

Indeed, by fitting the theoretical gaps to the experimental ones (Sec.IV.2) we have 
estimated the interaction parameters of the system. With these parameter values we 
have calculated other normal-state and superconducting properties. We have set up the 
formalism for the calculation of the particle-particle excitation spectrum in the RPA for 
the extended-Hubbard model in the CDW phase. In addition, we have demonstrated 
that a smooth interpolation between the large- and small-t/ theories is possible at T = 0, 
and that the particle-particle excitation spectrum we calculate at intermediate C/, evolves 
smoothly into the correct pseudospin-wave spectrum in the large-C/ limit. For typical 
intermediate values of f/, and V (i.e., when U < 2zt and V < U) our phase diagram in 
the T—x plane (Fig.IV.la), has the same qualitative features as the large-C/ phase diagram 
(Fig.III.l) and shares many features with experimental phase diagrams. Furthermore, the 
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charge-contrast in the CDW phase gets reduced and the superconducting coherence length 
enhanced over the values obtained in the large-f/ limit, in slightly better accord with the 
experiments. We have also argued that the processes that give rise to the attractive 
interaction, if electronic in origin, must become less and less effective as the system is 
doped towards the metallic side, so the strength of the attractive interaction ought to be 
reduced in order to obtain the experimentally observed superconducting Tc. 

It is worth stressing that our theory disfavours a purely phonon-mediated mechanism 
for the negative U (unless an alternative explanation can be found for the well-separated 
optical and transport gaps). A phonon mediated-attraction is necessarily retarded (cut off 
at characteristic phonon frequencies ~ 100 K); the minimum energy range over which the 
attraction must operate, in order to produce the two gaps, is about leV. In addition, the 
large binding required to produce the two-particle bound state at about 0.24eV above the 
chemical potential cannot come from phononic mechanism. It has been observed recently 
(Hinks, et al., |]14|, |l5l) that calculations based on realistic phonon these systems can, at 
least at the level of molecular-dynamics simulations, produce rather large superconduct- 
ing Tc, without requiring a very strong electron-phonon coupling {Xe-p around 0.6-0.8) 
or a large DOS at the Fermi level. The correct scenario could be that in the semicon- 
ducting phase of these systems, the attraction is mainly electronic in origin; however, 
this electronic component of the attraction becomes weak as the system is driven towards 
the metallic side (see discussion above), and the electron-phonon induced attraaction also 
becomes important in the superconducting phase. 

According to our theory the carriers responsible for charge transport in the semi- 
conducting phase of the bismuthate systems are charge-2e bosons (Cooperons). Hence 
experiments that can pick up the charges of carriers in the semiconducting phase would 
be of great value in verifying our theory. Unfortunately the determination of the charge 
of carriers in pure transport measurements is very difficult; they generally yield only the 



80 



ratio of the effective charge and effective mass. 



However, the following is worth looking into. Stanton and Wilkins ||104|| have observed 
that the noise spectrum of the high-field transport in a semiconductor is S{u!, E) = 
(ne*^ro)/(l + (t^To)^) x {A/2L)[kBTo/m* + (e*£'ro/m*)^], where E is the external electric 
field, To is the equilibrium lattice temperature, tq is the phenomenological relaxation time 
of the charge carrier of effective charge e* and mass m*, and A and L are respectively 
the area and length (in the direction of transport) of the sample. Now if the thermal- 
and electric-field dependent components of the noise spectrum could be separated, the 
above expression for S{uj, E) would allow us to calculate the charge of the carriers and 
their effective mass separately. To measure the noise spectrum, however, one must do 
a very controlled experiment at low temperatures on a pure, single- crystal sample of 
BaBiOs. A measurement of the phase change on passing through a magnetic field (Bohm- 
Aharanov effect) can also, in principle, detect the charge of the carrier, but such an 
experiment on mesoscopic samples of BaBiOs may be difficult to perform and hard to 
interpret. Even more likely possibilities for observing the charge-2e Cooperon may be in 
tunnelling measurements from a superconductor into the BaBiOs semiconductor, since 
pair tunnelling should now become viable. 

We believe that such experimental work that confirms or rules out this picture of the 
semiconducting phase of the barium bismuthates would be of great interest. A confirma- 
tion of this transport mechanism in the semiconducting phase of the barium bismuthates 
would confirm the electronic origin, as opposed to the phononic origin, for the negative 
U. 

We note that the phase diagram (Fig.IV.l), obtained for the negative-t/, extended 
Hubbard model in our mean- field theory does not contain any intermediate phase (i.e., 
one with both CDW and superconducting order), unlike the mean-field phase diagram of 
Micnas, et al ||80|. In the large-t/ case earlier work (Liu and Fisher and Matsuda and 
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Tsuneto |^) has shown clearly that, in mean- field theory, such a phase is possible only 
in the presence of longer range of interactions. 

Needless to say, in spite of the attractive features discussed above, the negative-t/ 
extended-Hubbard model, as applied to these bismuthates, has its limitations. The 
main shortcomings arise because we have neglected electron-phonon interactions,the long- 
range part of the Coulomb interactions, and disorder (charge disorder in the case of 
Bai^^.i^'a-BiOa). We comment on these shortcomings below. 

The neglect of Coulomb interactions is probably the main reason for one important 
discrepancy between the experimental phase diagram (Fig.I.2b) and the one obtained 
(Fig. IV. 1) from the model discussed above: The former has no two-phase region for 
X < 0.4, unlike the latter. Since the CDW and the superconducting phases have different 
average Bi-occupancies (or charges), long-range. Coulomb interactions strongly disfavour 
any bulk phase separation. Instead they might favour the formation, for small deviations 
from half-filling, of a new phase consisting of a dispersion of globules of superconducting 
material (which in the \aige-U limit and with small enough globule sizes can be thought 
of as discommensurations,citeCMV) inside a semiconducting CDW background. If the 
background compensating charge (because of potassium ions) were smeared uniformly, one 
might expect this dispersion to be a periodic array. However, in the actual system the 
potassium ions are distributed randomly, so the dispersion might well be random. Such 
a phase would be semiconducting above a new critical temperature (much lower than the 
critical temperature for the bulk superconductivity) corresponding to the destruction of 
the phase coherence between the different superconducting globules ( this temperature 
would be determined by the weak Josephson coupling between the globules). It is inter- 
esting to speculate that the semiconducting region, observed experimentally (Fig. III. 2b) 
away from half-filling, corresponds to such a dispersed phase. 

Another possibility that one ought to explore, both experimentally, and theoretically. 
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even in the nearest-neighbour model (IV. 2) and the mean-field framework, is whether the 
first-order coexistence region Fig.IV.lb gets replaced by an incommensurate phase, or a 
succession of commensurate phases. 

We must also account for the effects of disorder: In the potassium-doped case since 
K^"^ replaces Ba^"*", in addition to one electron being removed from the system (which 
we have considered in our model), a random Coulomb potential is also introduced (K is 
effectively a 1~ impurity). In the simple theory presented here, we have neglected this 
random potential. The effect of randomness is even more serious in lead-doped BaBiOs 
since the Pb atoms remove the negative-?/ centers of Bi randomly. Only very simple 
studies of the effects of this randomness on CDW ordering have been carried out so far 
(Jurczek and Rice and Jurczek |l83|). 



Since BaPbOs is a metal, it is clear that for a proper study of the entire range of lead 
doping, the one-band model (II. 2) is inadequate and the more complex model (II. 1) has 
to be used. Also the electron-phonon interactions are strong in these systems and has to 
be treated in a better way than has been done so far, especially because of the structural 
transitions in these systems. Clearly more detailed experimental and theoretical studies 
are called for to address these issues. 

It is worth noting in the end that there are some other theoretical attempts recently 
using negative-?/, extended Hubbard model to understand the bismuthates. Ohkawa 
||105|| has studied a negative-t/, extended-Hubbard model using 1/d expansion in the slave 



boson approach [106]. The leading order effect in l/rf produces a "kondo" like "charge 



fluctuation" regime discussed earlier in the context of an attractive impurity model |S| 



A mean-fleld analysis of this model produces, at a lower temperature, both CDW and 
superconducting phases as we have shown above rU|. Ohkawa argues that the normal 



state of barium potassium bismuthate represents a diamagnetic heavy electronic system 
with a large (~ 3 — 4 x lO'^K) charge-kondo temperature Tk- Superconductivity occurs 
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in a band of width 4Tk in the usual BCS route (via Cooper pairing), with the enhanced 
mass of the electrons leading to high Tc. A three band, negative-?/ extended Hubbard 
model has also been considered by Bala and Oles ||107|| recently leading to a phase diagram 
similar to ours (as described in Sec. IV). 
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Appendices 



Appendix A 

The Pseudospin Hamiltonian as the large-f/ hmit of the Hubbard Model 

There are various ways of transforming the Hubbard model to a spin model in the 
large-?/ limit. We use most straightforward route via second-order, degenerate, per- 
turbation theory, which is also physically transparent. Other methods use a canonical 



transformation 108, 109 



Consider the negative-t/ Hubbard model 



H = -t Y: c+c,. - Uj2in, - If. (Al) 

<ij>cr i 

In the large-?/ limit, the first term can be taken as a perturbation and one can construct 
an effective Hamiltonian describing the low-energy physics. 

The full Hilbert space consists of four states per site: the zero-occupancy, double- 
occupancy and single-occupancy (with up and down spins) states denoted by |0 >, |2 >, 
I |> and I |>, have energies —U, —U and respectively. The degenerate, lower energy 
states with double and zero occupancy are separated from the higher energy doublet of 
singly occupied states by the energy U. A chemical potential term lifts the degeneracy 
between the low energy doublet |0 > and |2 >. However, as long as f/ ^ yU, one can 
do a degenerate perturbation calculation involving the nearly degenerate lower energy 
manifold of unoccupied and doubly occupied states (see below). 

We now use degenerate perturbation theory to obtain the lowest-order (second order 
in t, as the first-order term vanishes) effective Hamiltonian in the subspace of states |0 > 
and 1 2 >, with the singly occupied states | t> and | J.> being virtual (intermediate) 
states. 
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It is clear that at this order, only the nearest-neighbour sites will be connected by the 
perturbation and so we choose a direct-product basis in the site representation: \iO;jO > 
, \i0]j2 >, \i2;j0 >, \i2]j2 > where |^0; jO > is the state with zero occupancy at both i 
and j sites, etc. The projected Hamiltonian in this subspace is 



H = 4^(|z0; jO X zO; jO| + \t2;j2 >< t2;j2\) + 4^(|z0; j2 X t2;j0\ + h.c). {A2) 
Now we define the operators 



St = \i2 >< iO\, 

S- = \i0><i2\, {A3) 

_ {\i2 >< i2\ - \iO >< iO|) 
Si = ^ , 

where \i2 > represents the state with 2 electrons at the i-th. site and Sj^ = {S^Y . 

These operators obey the commutation relations of spin operators and hence are called 
pseudospin operators. Note that 5*^^ = | imphes that the i-th site is a doubly occupied 
site while S'f = — | implies it is unoccupied. 

In terms of these operators the effective Hamiltonian becomes 

Hpseudospin J ^ ^ (Sj.Sj — ), (^4) 

<ij> 

where J = At^/U. 

This is an isotropic Heisenberg Hamiltonian. Deviations away from half filling, mod- 
elled by a chemical potential fi U, will break the pseudospin-rotation invariance 



and lift the degeneracy between the states |0 > and |2 >. So long as t/ ^ /i, the correc- 
tions to the results of degenerate perturbation theory because of e. they are 
clearly higher order in ^. The term VJ2<ij>^iTT-j, when projected onto the subspace of 
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|0 > and |2 >, takes the form 4:V J2<ij>{S- + ^){Sj + which also breaks the pseudospin 
rotation symmetry. 

Thus the final effective Hamiltonian in the projected subspace is 

Hpseud^spin ^JT. (S.-S, - |) + 41/ X: StS^^ - B Y,{2St + 1), (A5) 

<i3> <ij> i 

where B — /i + U — 2zV and z is the coordination number of the lattice. 
Appendix B 

The Relation Between Pseudomagnetic Field and Magnetization in the Dif- 
ferent Mean-Field Phases and the Phase Boundaries Across Them 

To obtain the pseudomagnetic field (and hence the chemical potential) in the different 
mean-field phases obtained by the mean- field theory of the pseudospin model (Sec. III. 2), 
we use the following thermodynamic identity 

-df/dB = n=<S^> + <S]> +1. {Bl) 

The subscripts i and j label sites on the two different sublattices A and B (or B and 
A). This equation, in conjunction with the self-consistency equations (Sec. III. 2) yield the 
pseudomagnetic field in in terms of the magnetisation m. 

B.l Pseudomagnetic Field 

Superconducting Phase 

In the superconducting phase < S% >—< S% >=< > and | < 5"^ > 
Hence Ra = Rb = E (See Section III). Thus equation (III. 2) gives 



<S%> 
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where m — n — 1, is the magnetization, Jq — zJ and Kq — zK. 
Prom the self-consistency equation (Sec. Ill) 

< Sl^B) >= < S%^A) > tanh(/3i?B(A)). 

For the superconducting phase 

^„-' = m 

By combining equations (B2) and (B3), we get 

B^{Ko + Jo)m/2 (BA) 

in the superconducting phases. 
Nonordered Phase 

Similarly for the nonordered phase (< S'^(b) >= and < -S"^ >=< >=< > ) 
whence 

k 

B ^ 1 + p-Hanh-^ m + ^m. {B5) 

Chcirge-Density-Wave Phcise 

Here < -S"^ >=< S% >— 0. So, following the procedure outhned above we get 

m = l/2[tanh/?(S - Ko/2{m - -q)) + tanh/3(B - KQ/2{m + 77))], (56) 
where 77 =< S\> - < S%> 



B.2 Phase Boundaries 



Charge-Density- Wave — Nonordered 

The CDW-nonordered phase boundary is obtained by noting that the staggered mag- 
netisation ?7=<5'^> — <5'^> goes to zero continuously as this boundary is approached 
from the CDW side. By expanding the self-consistency equations about this boundary 
(?7 ^ ) we get 

kBTcN = ^{l-m'). (57) 

where Tqn is the transition temperature. 

Superconductor-Nonordered 

Here the order parameter < 5*^' >^ as the phase boundary is approached from the 
low-temperature side. An expansion of the self-consistency equations yields 

B = Kom/2 + Jotanhp{B - Kom), (58) 
which along with Eq.(B4) gives 

TsN = Jom/ikB log[(l + m)/(l - m)]). (59) 

Appendix C 

The Equation of Motion for the Green Function in the RPA 

To obtain the excitation spectra in different mean-field phases, we used a Random- 
Phase- Approximation technique ^ for the closure of the equations of motion for the 
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Green function (Sec. III). We use the Hubbard operators | |100[ | P^^ = \ia >< In 
the models we discuss a and /? can take two values (+ and -) corresponding to the two 
mean-field eigenstates |+ >= cos|| t> +sin|| J,> and |— >= cos|| l> — sin|| t>. 
The indices i label sites on a bipartite lattice; with nearest-neighbour interactions, the 
only relevant indices are the two sublattice labels (A and B). The commutation relation 
between these Hubbard operators is 



{CI) 



where [R,S]_ represents a commutator between R and S. In terms of these projec- 
tion operators the double-time, retarded Green function G^'^^^it — t') (also denoted by 



G7p,,sit - 1') = -^oit - 1') < [Pm^ m')]- > ■ (C2) 

This Green function gives the probability amplitude for the process in which, at time 
t', the n*'* site makes a transition from the state \n6 > to the state {wy > followed by 
another transition, at site m and a later time t, from the state \m(3 > to the state \ma >. 
The Fourier transform of this Green function in time is 

GX^M = ^ r rftGX,(t-t')exp.c.(t-t'). (C3) 

ZTT J — oo 

Correlation functions of two projection operators are expressed in terms of the Green 
function in the usual way: 

POO 

J —oo 

where fsi^) is the Bose occupation factor {expu/ksT — 1)^^ 
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Since any local operator can be represented in terms of the projection operators as 
A = Sa/3(A)o!/3-Pa/3' psBudospin Hamiltonian can be written as 



where 



and 



ul^ = B[2{Sn^p + l] 



An obvious symmetry of W^/^^^g, which we use, is W^^^i^p^, — W^^t^p^,. The equation 
of motion is, therefore 



idGt^',pp'{t - t')/dt = ^S{t - t') < [PL',Pk'^- > + « [PL', HUPU » ■ m 

If we substitute Eq. (C5) in Eq. (C7), rearrange dummy indices , and use Eq.(Cl) 
and the symmetry of W^^i^^p,, we obtain the following equation of motion for the Fourier 
transform (i.e., in q, a;-space) of the Green function. 



+ E « PU'P'APk' » ■ iC8) 

k,iJ,vv' 

In terms of the Hubbard operators, the RPA is defined by the decoupling in which a 
three-operator Green function is factorised into products of two-operator Green functions 
and expectation values of the remaining operator, as shown below: 
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<< Paa'P^-y'WJjji' »oj— Saa' < Paa' > ^77' ,/3/3' ('^ ) + ^77' < -^77' > ^aa' ,/3/3' ) 

In our model (C5), the last term on the right hand side of Eq.(C9), arising out of the 
commutation of P^/ and Pl^^i , vanishes as there is no interaction connecting the sites on 
the same sublattice (the delta- function 5im requires this but W^^/^^^z = for i = j). In a 
model with interactions connecting the same sublattices (e.g., a model with further than 
nearest-neighbour interactions) this term has to be considered. 

By using the RPA factorisation in Eq.(C8) we obtain 



k,tJ,i/ 



where 



e'a = + E W^a,77(q = 0) < ^77' >' (^H) 

1,1 



and 



Since the site indices can have 2 values (denoting the 2 sublattices), and since a and 
a' can have only two values, the most general form of the above equations of motion here 
is a 4 X 4 matrix equation given by 
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/ 


UJ - 2Rb 





^^+q 






( ^+-13/3' 


:q,^) \ 







uj + 2Rb 


-q 


-n^q 




^-+13(3' 










UJ - 2Ra 







^+-1313' 




V 


q 


q 





iJ + 2RA ) 




I pBB 


:q,^) / 



27r 



-L 



B 



( \ 



5+/3'5/3- 



{CU) 



where 13 [3' take values H — or — h and 



n!tq= -[±J(q)(l + cos^^cos0B) -ir(q)sin^^sin^ij]LY_. (C15) 

The sublattice indices can be dropped in the superconducting and nonordered phases 
(i.e., there is no broken lattice-translational symmetry). This matrix equation (C14) 
reduces to a 2x2 matrix equation in the CDW and the superconducting phases and to 
a scalar equation for the nonordered metallic phase because of the symmetry of these 
phases. 



Appendix D 

Scaling forms for Thermodynamic Functions near Multicritical Points 

In the phase diagrams of Figs. III. 1 and III. 2 the following multicritical points occur: 
bicritical points, tricritical points, tetracritical points and critical end points. For the sake 
of completeness we give the scaling form [ |1 1 1| for the singular part of the free energy in 



the vicinities of these points. 
Bicritical and Tetracritical Points 



Bicritical and tetracritical points are studied easily using the same model, namely, the 
Landau- Ginzburg free-energy functional [^, ^ p.l2| . 
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i=l 



where ip is a three-component, vector order parameter with components '0,, z = 1 . . . 3. A 
bicritical point is obtained if f < and a tetracritical point i{v>0. The temperature-hke 
scahng variable t depends hnearly on the deviation of ro from its value at the multicritical 
point. 

The singular part of the free energy JF^ assumes the following scaling form in the 
vicinity of a bicritical point which occurs at the temperature Tf,: 

where the mean-field values of the exponents are = 0, 0^ = 1 and 0^, = 0. The specific- 
heat exponent at is related to the correlation length exponent via the hyperscaling 
relation dub = 2 — at, where d is the dimension ( the mean-field value of 1/5 = 1/2 ). To 
leading order in e = 4 — ci the bicritical exponents are 

= 1/2(1 + 5e/22), 

(j)g = l + 3e/22, {D3) 
(j), = -e/22. 

Note that, even though v is irrelevant, it should be displayed explicitly since it is 
a dangerous irrelevant variable. More details regarding this point and a comprehensive 
discussion of the (similar) scaling in the vicinity of tetracritical points can be found 
elsewhere |9|, [m], ITTl . 



For the anisotropic antiferromagnets, that we are interested in, the Landau- Ginzburg 
free energy is an expansion in powers of M and M^, the magnetization and the staggered 
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magnetization. However, in the vicinity of the bicritical point, such a free energy can 
always be put in the form (CI) given above. 



Tricritical Point 

In the vicinity of a tricritical point the Landau- Ginzburg free-energy functional can 
be written as m, |lT|, [USj ITTGI, [TT7 



where ip is a scalar order parameter. The tricritical point occurs at = 0, 02 = 04 = 0, 
and Og > ( a2 depends linearly on temperature ). The scaling form for the singular part 
of the free energy J-'s has the following form in the vicinity of the tricritical point Tf. 



g JJ 

where, in mean-field theory, the exponents are at = 1/2, (pt = 1/2, At = 5/4, and the 
scaling field g is defined to be proportional to a^. Since the upper-critical dimension 
for tricritical behaviour is 3, corrections to mean-field theory are weak logarithmic ones 
\TmKT3\. 



Critical End Point 

The simplest Landau-Ginzburg theory which yields a critical end point is defined by 
the functional Im], |n8|, |1T|, ^ 



with V < 0,w > 0, and a scalar order parameter. (For the critical end point 
shown in Fig. III. 2 a more complicated Landau-Ginzburg functional would be required, 
but Eq.(D6) suffices for the purpose of illustrating the essential features associated with 
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this point.) If ip is uniform and we use mean-field theory (i.e., we minimise J^lg with 

2 

respect to ip)^ then a critical end point is encountered at m = ^ and r = 0. As this point 
is approached from the disordered phase, thermodynamic functions, such as the order- 
parameter susceptibility Xi diverge as they do along the critical line which terminates at 
the critical end point (Fig.III.2); i.e., at the mean-field level X'S'j~r~^ as r ^ (as usual r 
depends linearly on the deviation of the temperature from the temperature at which the 
critical end point occurs). However, the order parameter jumps from zero to (|f |/2u;)^/^ 
as we go from the disordered phase to the ordered phase through the critical end point. 
This behavior does not depend on the mean-field approximation. Even renormalization- 
group calculations [ |118| , |119| , |120|| yield the same qualitative behavior: Thermodynamic 
functions like x diverge with the same exponent at a critical end point as they do along 
the critical line which terminates at this point and the order parameter jumps at the 
critical end point. 

Appendix E 

The Current Operator in the Pseudospin Model 

The effective Hamiltonian in the pseudospin model is (see Appendix A) 

^ = E st.s^, +4-jT. ^^ts- + B.C.) + wy: s!.s] 

<ij> <ij> <ij> 

-5^(25^ + 1). (£^1) 

i 

To couple this to an electromagnetic field we use the Peierls substitution t 
texp—i // A.cil where A is the vector potential {Ti = c = 1). Such a substitution will 
obviously leave the first term in H unaffected and the only contribution to the current 
comes from the transverse part of the pseudospin-pseudospin interaction term. 

By expanding H to linear order in A, and noting that a minimal coupling to the 
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current is of the form j. A, we can easily identify the current operator for the above model 
as 

j = -i{2e) Y: MS^.Si - H.C.), {E2) 

<ij> 

lip 

where a is the nearest-neighbour vector and J — 4'-^. 

There are alternative ways of deriving the current operator (of course with identical 
results) for the projected Hamiltonian . For example, we can construct the dipole operator 
P for this model and then use the relation j = — i[P, H]_ to get the current operator. 

Appendix F 

The Ginzburg-Landau Equation and the Superconducting Coherence Length 
for the Pseudospin Model 

The Ginzburg-Landau free-energy functional for the pseudospin model in terms of the 
order parameters < S^' >— V'i and 2 < 5"? >= (n — 1) = m is: 

J^[^] = -J/N {4^j + H.C.) + Xom^ 

<ij> 

+{PN)-' YiiSi - 1/2) log(5, - 1/2) + {Si + 1/2) log(5, + 1/2)], (Fl) 

i 

where N is the number of sites, and 

Sf^m' + my^. (F2) 

It is then easy to show that an expansion to order ip'^ gives the phase boundary between 
the superconducting and nonordered phases. To this order the free energy per site is 



= -J/N Y ii^li^j + H.C.) + {(3N)-' Y mV{mlogn/{l - m)). 

<ij> i 



(F3) 
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which of course gives the same expression for Tcn as Eq.(B8) 

Note that the couphng to the electromagnetic field is expressed in terms of the following 

gauge transformation on the bond (ij), 

ipiil^j ijjiijjj exp -i(2e) J^^ A.dl. (F4) 

where A is the vector potential. Note that the current operator (Appendix E) follows 
from a functional derivative of J-' with respect to A 

To calculate the coherence length, we take the magnetic field H to be in the z direction, 
so that, in the Landau gauge, the vector potential is A = (0, —Hx, 0). Then 

where 

= 2eHl^a''. (F6) 

Here a is the lattice constant and Ix is the x-coordinate of the l-th site; 

Close to the transition, by expanding the Ginzburg- Landau free energy and retaining 
terms upto order of -0^, we obtain 



J^ii^] = -J/N ^ (V^lV^,- + H.C.) + T/Tsn{Q) E 

<ij> i 

where T5iv(H) is the transition temperature in a magnetic field H. 

An expansion of the first term on the right hand side gives the term. By 

equating the coefficient of this term to (where ^ is the superconducting coherence 
length) we obtain 

e = a'/r, (F8) 
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where r = 1 — T/Tsiv(0). ( Equivalent ly, one can write down the equation for ipi and 
compare it with the eigenvalue equation of a quantum oscillator.) 

Hence, as wc have already mentioned in the text, the coherence length well away from 
the transition is about one a lattice constant. This result is meaningful from a physical 
point of view as, in the large- C/ theory, there is no length scale in the problem other than 
the lattice constant. 
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Note 1. Some other correlation functions are gapless in the superconducting phase 
reflecting merely the XY-model nature of the ordering; however, they do not correspond 
to any physically measurable function. 



Note 2. In an earlier treatment |PD| it was claimed that a mixed phase with both 
CDW and superconducting order could be obtained within mean-field theory for model 
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by working with fixed < > (the analogue of the canonical ensemble). This is not 
correct. The phase diagrams obtained from both fixed--B (i.e., grand canonical ensemble) 
and fixed-< > (i.e., canonical ensemble) calculations have to be the same by the 
equivalence of the two ensembles. The first-order boundary of Fig. III. la (obtained by 
controlling B) appears as the two-phase coexistence region of Fig.III.lb (obtained by 
controlling < >); and the discontinuity of < S*^ > across the the first-order boundary 
of Fig.III.lb yields directly the boundaries of the coexistence region of Fig.III.la. The 
only difference between the two ensembles is the following: At coexistence and with fixed 
< 5^ > , the system phase separates into two phases-one superconducting and the other 
a semiconducting, CDW phase-with one interface between the two. (The fractions of 
the two phases follow from the lever rule.) However, with fixed B on the first-order 
boundary), either one or the other phase fills the whole system (since < > is not 
fixed) and there is no interface; formally one or the other can be picked out by allowing 
B to approach the first-order boundary from either the superconducting or the CDW 
side. As we have discussed in the text, to obtain a bona fide mixed phase with both 
super CO ducting and charge-density-wave order, next-nearest-neighbour interactions must 
be added to the pseudospin model. This has also been pointed out recently by . . . 

Note 3. The bose condensation of the hard-core bosons is observed from the behavior 
of the chemical potential as one approaches the nonordered-singlet superconductor phase 
boundary from the high-temperature side. The chemical potential loses explicit depen- 
dence on temperature on entering the superconducting phase from the high-temperature 
disordered phase signifying the onset of a bose condensation of the hard-core bosons. 
The pseudomagnetic field in the singlet superconducting phase is given by [Eqn (B4)] 
B = {Kq + Jo) < 5^ > and is clearly independent of temperature. 

Note 4. From the width (~ 6 eV) of the antibonding band of the LAPW, band- 
structure calculation [^, The value of the hopping matrix element is obtained from 
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the experimental values of the (renormalized) DOS at the Fermilevel when fitted to 
a tight binding band with square density of states, it is less (~ 0.20 eV) than the value 
(0.3 eV) obtained from the band-structure calculation. 



Note 5. In earlier work [ 101 | the only a singlet superconducting order parameter b^. 
In a negative-?/, extended-Hubbard model-context, choice is naive as we have shown here; 
both b^ and be are nonzero in the superconducting phase. 

Note 6: This symmetry is quite easily understood: The negative-t/ Hubbard model 
can be transformed into the corresponding positive-?/ model, as described earlier. (If the 
chemical potential of the negative-?/ model is non zero (i.e., the system is away from half 
filling) this transformation leads to a positive-f/, Hubbard model with a magnetic field 
coupling to the spin.) The positive-f/ Hubbard model has a global SU(2) spin rotation 
symmetry that mixes the up and down spins (Affleck, et al. (1988)), whence the x{y) 
ordered and the z ordered antiferromagnetic (SDW) states are degenerate. In the negative- 



U context, this symmetry is a pseudospin rotation symmetry which means that the 
superconducting state is degenerate with the CDW state at all filling. [However, as the 
negative-?/ model away from half filling maps onto a positive-?/ Hubbard model with a 
magnetic field, this symmetry between the x{y) and the z-ordered states is lifted, i.e., the 
corresponding symmetry in the negative-?/ model, between CDW and superconducting 
states is destroyed away from half filling. 



Ill 



Table I.l 



Normal State Properties of BaPbi_j;Bij;03 



Metallic 


Semiconducting 


X 


0.0 


0.12 


0.2 


1.0 


p(i7cm) X 10~^ 


2.10 


4.15 


5.40 




n(cm~^) 


1.42 X 10^0 


16.0 X 10^0 


20.0 X 10^0 


Prefactor 110 x 10^° 
(activated with gap 
0.24 eV) 


/x(cm^v~^s~^) 


210.0 


9.40 


5.80 




Amf(^) 


220.0 


22.0 


15.0 




m* /rrie 


0.16 


0.32 


0.23 





Table 1.2 



Superconducting Properties of BaPbi-^^Oa and Bai-j^K^^BiOa 





Bao.6Ko.4Bi03 


BaPbo.75Bio.25O3 


dHc^/dT 


-5±0.5KOe/°K 


-5.3KOe/°K 


dHc/dT \rp^ 


-75.0 ± 8.0 Oe/°K 


-50±5Oe/°K 


-Xgl 


45.01 


80.01 




70.0A 


80.01 


AC/T, 


2.1 ±0.4mJ/moleK^ 


0.94 mJ/moleK^ 


iV(0) 


0.42 states/eV 


0.42 states/eV 




< 0.64states/eV 


< 0.64states/eV 


^cl-ph 


= 0.5 , (but less than 


= 0.7, (but less than 




0.8) 


0.8) 


7(= AC/1.437;) 


^ 1.5mJ/moleK^ 


0.8 ± 0.2mJ/moleK^ 


-^penetration 


^ 4000A 


5000 - 6OOOA 
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Table II. 1 



U3+ = 19.7 U3+ = 18.8 

Bi U4+ = 19.7 Sb U4+ = 11.9 
U5+ = 32.3 U5+ = 52.0 



U2+ = 16.9 
Pb U3+ = 10.4 
U4+ = 25.5 
Ui+ = 14.3 
Tl U2+ = 9.4 
U3+ = 20.9 
U3+ - 13.3 
Nb U4+ = 11.7 
U5+ = 53.0 



U2+ = 15.9 
Sn U3+ = 9.1 
U4+ = 31.6 
Ui+ = 13.1 
Sn U2+ = 9.2 
U3+ = 26.4 
U3+ = 18.7 
V U4+ = 17.0 
U5+ = 64.0 



U2+ = 18.2 
Ge U3+ = 9.5 
U4+ = 48.7 
Ui+ = 14.6 
Ga U2+ = 9.2 
U3+ = 26.4 
U2+ = 13.9 
Ti U3+ = 15.8 
U4+ = 56.6 



Table II.2 



Correspondence between Pseudospin model & Extended Hubbard model 



Pseudospin model 


Extended Hubbard model 


Up spin 


Vacancy 


Down spin 


Double occupancy 


Magnetisation m 


Deviation from | filling 


Magnetic field B 


Chemical Potential fi 


Antiferromagnetic (AFM) order 


Bose Crystal of Cooperons (CDW 
phase) 


Spin flopped phase (i.e. AFM 
with nonzero XY component of 

spins) 


Bose condensate of Cooperons (S- 
wave superconductivity) 


Non-ordered phase (with partial 
alignment for B ^ 


Bose liquid of the Cooperons 
(normal metal) 
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Figure Captions 

Fig. I. la The underlying cubic perovskite structure of BaBiOa with oxygen atoms at the 
corners of the octahedra centered around Bi-atoms. Bi-atoms occupy the vertices 
of a simple-cubic lattice; Ba atoms are at the body centres of this cubic lattice. 

Fig. I. lb Possible distortions of this underlying cubic perovskite structure that lead to 
the various lower-symmetry structures found in Bai-a^K^-BiOs and BaPbi-ajBi^-Os. 

Fig. 1. 2a Schematic phase diagram of the BaPbi_a;Bia;03 system in the T — x plane (after 
Sleight im). The orthorhombic(0//) and tetragonal (T) phases are metallic at high 
temperatures and superconducting at low temperatures. The orthorhombic (O//) 
and monoclinic(M) phases are semiconducting. The precise nature of the metal- 
semiconductor is not clear, but here it is shown as a first-order transition with the 
hatched region indicating a region of two-phase coexistence. All the four phases 
undergo high-temperature structural transitions to a cubic phase. 

Fig. 1. 2b Schematic phase diagram of the Bai-^K^BiOa system in the T — x plane after 
Pei, et al. [§, Q The entire cubic phase is metallic and the hatched, low temperature 
region is a superconducting phase. The monoclinic(M), orthorhombic (O), and 
rhombohedral (R) regions are semiconducting. 

Fig. 1. 3a Temperature dependence of the resistivity of BaPb1_2.Bia.O3 on a log-log scale 
in the metallic regime and across the metal-semiconductor transition regime. Note 
that the temperature-coefficient of resistivity is negative at x = 0.3, but resistivity 
goes to zero at low temperatures. (Afte Uchida et al. p[) 

Fig. 1. 3b The temperature dependence of the electrical resistivity of BaPbi-^Bi^Os in 
the semiconducting regime (after Uchida et al. 0). 
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Fig. 1. 3c Temperature dependence of the electrical resistivity of Bai-^K^BiOa (for x = 
0.45) thin films for two different samples i and j prepared under slightly different 
annealing conditions (see Sato, et al. [^). 



Fig. 1. 4 Reflectivity spectra in the metallic regime < x < 0.2 of BaPbi-ajBi^jOs alloys. 
The solid curves are results of a Drude fit. Note the deviation of the Drude fit a 
X = 0.2 (after Uchida et al. 0). 

Fig. 1. 5 Variation of the optical gap Eg (i.e., the position of the peak in the optical 
absorption spectrum) and the transport gap Ea (i.e., transport activation energy) 
in eV, of BaPbi_a;Bia;03 with concentration x, in the semiconducting phase. The 
optical gap does not vanish at the semiconductor-metal transition. Two different 
scales have been used to show the wide separation of the energies of the two gaps 
(after Uchida, et al. g). 

Fig. 1. 6a Optical conductivity spectrum (in a linear scale) in the semiconducting phase 
of BaPbi_xBi2.03 alloys. Note the gradual shift of the conductivity peak towards 
high frequency with increasing x (after Uchida, et al. 0). 

Fig. 1. 6b Optical conductivity of Bai-^^K^-BiOs with x = 0.0 (solid), x = 0.1 (dot), 
X = 0.2 (dash), x = 0.33 (dot dash) and x = 0.4 (solid) from Blanton, et al. . 



Fig. 1. 7 (a) The composition dependence of the magnetic susceptibility at 150K in 
BaPbi_2.Bi^03. The solid line is the diamagnetic contribution to susceptibility 
from the ionic core electrons (after Uchida, et al. [^). In (b) is shown the normal 
state magnetic susceptibility of both lead and potassium doped systems as a func- 
tion of temperature. The diamagnetic core contributions for the lead and potassium 
doped systems are, respectively, 8.6x10^^ emu/mole and 7.8x10"^ emu/mole (after 
Batlogg, et al. pi). 
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Fig. 1. 8a The temperature dependence of the resistivity of a single crystal thin-film sam- 
ple of BaPbo.7Bio.3O3 at different levels of oxygen deficiency. The labels 5 to 1 are 
in the increasing order of oxygen deficiency (after Enomoto et al. [|63|1 ). 

Fig. 1. 8b The temperature dependence of the carrier concentration for a BaPbo.7Bio.3O3 
single-crystal thin filsm. The labels 5 to 1 are in order of increasing Oxyen deficiency 



[after Enomoto et al. |63 



Fig. 1. 9 Temperature dependence of the upper critical field (in Tesla): (a) For 
BaPbi_2.Bi^03 at different filling x. S and P denote, respectively, single crystal 
and polycrystal samples. Doping x and the superconducting Tc of the samples are 
given with each curve (after Kitazawa, et al. [0). (b) For Ba1_2.K2.BiO3 thin films 
at a; = 0.45 (after Sato, et al. jiSj). 



Fig. 1. 10 A spherically symmetric bismuth 6s orbital surrounded by the oxygen p orbitals 
in the BiOe octahedron. Only the a orbitals of each of the six oxygen atoms are 
shown. 

Fig. 1. 11 LAPW band-structure results for BaPbi_2;Bi2.03 alloys. The effect of Pb doping 
is shown from (c) through (a); a and a* denote the bonding and antibonding a bands 
respectively. The antibonding band that crosses the Fermi level Ep comes from the 
Bi/Pb6s-02pCr bonding. Lead doping adjusts the filling of the a* subband and also 
introduces a noticeable chemical shift in the position of the 6s band (open circles) 
(after Mattheiss 



Fig. 1. 12 Comparison of schematic bands for the high- and low-T^ oxide superconductors, 
(a) The Fermi level lies in the antobonding band of predominantly p character for the 
high-Tc materials, (b) For the low materials the Fermi level is in the degenerate 
77* manifold. The quantities ep, e^, and e, are the energies of the hgand p-level and 
cation d and s levels, respectively (after Mattheiss 
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Fig. 1. 13 Band-structure results for the lead-doped system with a tight-binding fit. (a) 
LAPW energy bands (solid lines) - for cubic BaPbo.7Bio.3O3 compared with those 
derived from simple tight-binding models involving three (dotted lines) and five 
(dashed lines) tight-binding parameters respectively |Q. Note the ten bands are the 



high-lying Ba-derived bands, (b) Sketch of the electron-like F-centered simple cubic 
Fermi surface predicted by the three-parameter tight-binding model for BaBi03. 
Note that the electron surface nests with an identical hole surface centered at the 
Brillonin-zone corner (R) (after Mattheiss and Hamman 0). 

Fig. 1. 14 Model tight-binding conduction bands for cubic (dashed line) and distorted 
(solid line) BaPbi_2.Bia;03 alloys. The downward movement of the Fermi level as a 
function of doping is also shown. Note that at half filling the gap is at the Fermi 
surface (after Mattheiss p5|). 



Fig. 1. 15 (a) LAPW band-structure results for the Bai_2;K2;Bi03 system. The wide bands 
of Bi6s-02p are bounded above and below by the a* and a bands. The Fermi level 
lies in the a* band. Note the similarity between the band structures of the K-doped 



and Pb-doped systems (after Mattheiss and Hamman 0, ^). (b) The full potential 
LMTO band structure near fermi level for 2BaBi03. Solid line is for t = 0° and 
b = OA, dashed line is for t = 10° and b = OA, dot-dashed line is for t = 10° and 
b = 0.09A. The brillouin zone has been folded (after Liechtenstein, et al. p6|). 

Fig. 1. 16 Tc versus the Sommerfeld constant 7 for various superconductors, bismuthates 
and cuprates have the highest Tc's on an absolute scale (and especially given their 
low DOS at Ep). The extraordinary character of these two types of systems relative 
to all other superconductors, including oxides, is highlighted in the above plot. The 
dashed line demarcates low-Tc superconductors from the high-Tc superconductors 
(after Batlogg 0). 
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Fig. II. 1 Tc as a function of [/e///2V calculated in the Hartree approximation for the 
BEG model (lower curve). Solid lines denote second order and dashed lines denote 
first-order transitions and the dot, a tricritical point. In the Hartree calculations Uo 
and e'^/aa were kept fixed at the values 2 and 1, respectively, and g'^/C is varied 
(after Rice and Sneddon |]78| ). 



Fig. II. 2 Tc as a function of z, the average coordination number in the BaPbi-^^^Bi^Os 
system for three different values of f/e///2V. is calculated in the Bethe-Peierls 
approximation. The inset on the top right corner shows the effect of a nonzero 
band with {tij ^ 0) in this phase diagram, on the region near zero temperature, as 
described in the text, for Ueff/2V > 0. In this diagram M is a metallic phase and 
S stands for BCS superconductor (after Sneddon and Rice |7S|). 

Fig.II.3a The DOS of the Bi A-site (sohd hue) and Bi B-site (dashed line) at 20% Pb 
doping using eo = 2eV and A = 0.8. The Bi subbands approach each other because 
of a level repulsion between the Bi and Pb single bands. The total DOS is shown 
in the inset at the top right corner. The different level shifts of the Bi subbands 
produce a pseudogap around the Fermi energy. In (b) is shown the case of 60% Pb 
doping using the same parameter values as in (a). The Bi subbands continuously 
approach each other with increasing Pb doping. The states around the Fermi energy 
are still almost composed of Bi sites (after Jurczek and Rice ||79|| ) . 

Fig. II. 4a The boundary between semiconducting and metallic phases in the e'-^-' — x 
plane for several choices of U . The symbols are explained in the text. In (b) is 
shown the DOS obtained in the CPA at two different values of U (after Yoshioka 



and Fukuyama |89|| ). 



Fig. II. 5 The superconducting transition temperature as a function of filling x for different 
values of U with e'-^-' kept fixed. Notte that the phase diagram is symmetric with 
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respect to x — > 1 — x (after Yoshioka and Fukuyama 

Fig. III. la A typical phase diagram of the pseudospin model in the pseudomagnetic field 
(B) and temperature (T) plane. The full lines are second-order phase boundaries 
and the dashed lines indicate first-order phase boundaries. The two continuous lines 
and the first-order line meet at a bicritical point (the full circle). The first-order 
line separates the charge-density-wave (CDW) phase (antiferromagnet in the spin 
language) from the singlet superconducting (SS) (spin flopped in the spin language) 
phase. The two continuous lines separate the non-ordered metallic (NOM) phase 
(paramagnet in the spin language, with a net moment because of the aligning field) 
from the CDW and SS phases. The arrows show the arrangement of the spins on 
nearest-neighbour sites. The parameters used are Jo = 1.0 and Ko = 1.301. 

Fig. III. lb The phase diagram of Fig.III.la in the magnetisation (m) and temperature 
(T) plane. The full lines are second-order phase boundaries and the dashed lines 
indicate first-order phase boundaries. The hatched region is the two-phase coexis- 
tence region (SS and CDW) corresponding to the first-order phase boundary between 
CDW and SS phases in Fig.III.la. 

Fig. III. 2 The mean-field phase diagram when a next nearest neighbour interaction is 
included in model II. 5 (in the large-U limit). The phase denoted by I is the so 
called intermediate phase (after K.-S. Liu and M.E. Fisher, J. Low Temp. Phys. 
10, 655-683 (1973)). 

Fig. III. 3 The RPA excitation spectra of the pseudospin model shown along different 
symmetry directions of the Brillouin zone of a cubic lattice (see inset), with Jq = 
1.0 and Ko = 1.301. The lines marked CDW(+) and CDW(-) correspond to the 
upper and lower branches of the spectrum in the CDW phase (for temperature 
T = 0.25, pseudomagnetic field B = 0.25, magnetisation m = 0.4 and staggered 
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magnetisation rj = 0.96). The dashed hne shows the spectrum for the SS phase 
(at T = 0.25, B = 0.70, m = 0.61 and Oa ^ Ob ^ 50.55°, where 9a & Ob are 
defined in the text). The dash-dotted hne shows the spectrum in the NOM phase 
(at T = 0.80, -B = 0.40 and m = 0.27). The spectrum in the SS phase is gapless 
whereas the spectrum in the CDW phase has the anisotropy gap discussed in the 
text. 

Fig.III.4 The optical conductivity (in arbitrary units) in the large-C/ hmit using RPA 
plotted as a function of frequency. The parameter values taken are zt = 1 eV, U = 2 
eV, and zV = 0.05 eV. The same phenomenological broadening (F = 0.09) has been 
used for both the channels. Note that the two-particle absorption peak is a factor 
of about 9 down compared to the one-particle peak. 

Fig.IV.l(a) The mean- field phase diagrams of the negative-U, extended Hubbard model 
in filling-temperature plane, for two different values of the interaction parameters. 
The solid lines and the hatched region have the same meaning as in Fig. III. 1. The 
values of Coulomb interactions and zt used to obtain this phase diagram are shown 
in the figure. 

Fig. IV. lb The mean-field phase diagram of the negative-U, extended Hubbard model 
in the chemical potential-temperature plane. The solid and the dashed lines denote 
continuous and first-order transition, respectively as in Fig.III.la. Other symbols 
also have the same meanings as in Fig.III.l. 

Fig. IV. 2 (a) The single-particle excitation spectrum in the mean- field CDW phase where 
Ac is the CDW gap parameter. In (b) are shown, schematically, the two-particle 
(hole) bound states that form inside the gap of the two-particle spectrum. The 
energy of these bound states is the transport gap {Ea) in our theory. 

Fig.IV.3 The diagrammatic representation of the RPA series for the pair susceptibility 
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in the mean-field CDW state. Double lines denote mean-field single-particle Green 
function in the CDW state and the wavy line denotes the interaction ^'jj;i'j'(k,k')- 

Fig. IV. 4 A schematic representation of the Cooper instability equation as described in 
the text. 

Fig.IV.5a The bare (empty circles) and the screened (empty triangles) onsite Coulomb 
interactions as functions of (half) the band width. The screening mechanism is 
explained in the text (see also Fig. IV. 7). 

Fig. IV. 5b The extended Coulomb term V {x2z) as a function of (half) the band width 
W{— zt). The condition that we have used to obtain this graph is that the bound 
state energy should be equal to the transport gap (see the text). Note that in the 
bandwidth — limit the value of 2zV approaches 0.24eV, the result obtainted in 
the large- [/ limit. 

Fig. IV. 6 The Coulomb interactions U and zV as a functions of the bound-state energy. 
This graph is obtained with the condition that the CDW gap parameter is half the 
optical gap of BaBiOs (i.e., the CDW gap is 1 eV). Note that V goes to zero as the 
bound-state energy approaches zero. 

Fig. IV. 7(a) A diagrammatic representation of the screening (in RPA) of the particle- 
particle interaction. Note that the particle-particle interaction is screened by the 
particle-hole excitations. 

Fig.IV.7b Same as (a) for the particle-hole interaction. The particle-hole interaction 
is screened by the particle-particle processes. Note that our calculation is not self 
consistent (see text). 

Fig. IV. 8 (a) The superconducting transition temperature with and without the screening 
described in the text, (b) The CDW transition temeperature with and without the 
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screening. 

Fig. IV. 9 The supcrconductor-nonordered and tiic CDW -nonordered piiase boundary 
evaluated with the same parameters as in IV. lb. To show that screening of V 
does not much affect the superconductor-nonordered phase boundary at the range 
of parameter values we are interested in, we draw this phase boundary for V = 
also (dashed line). Note that this boundary is quite close to the boundary with 
2zV — 0.13eV (the solid and the dashed hne both represent continuous phase 
transitions.) 
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